Entanglement subspaces, trial wavef unctions, and special Hamiltonians 
in the fractional quantum Hall effect 
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We consider systems of trial wavefunctions in the fractional quantum Hall effect that can be 
obtained in various ways. In one way, the functions are obtained by analyzing the entanglement 
of the ground state wavefunction, partitioned into two parts. In another, functions are defined by 
the way in which they vanish as several coordinates approach the same value, or by a projection 
operator Hamiltonian that enforces those conditions. In a third way, the functions are given by 
conformal blocks from a conformal field theory (CFT). These different sets are closely related. The 
use of CFT methods permits an algebraic formulation to be given for all of them. In some cases, 
we can prove that all of these spaces are the same (a finite-size bulk-edge correspondence), thus 
answering several questions and conjectures. We can also use the analysis of functions to produce 
a projection operator Hamiltonian that produces them as zero-energy states. For a model related 
to the M = 1 superconformal algebra, the corresponding Hamiltonian imposes vanishing properties 
involving only three particles; for this we determine all the wavefunctions explicitly. We do the same 
for a sequence of models involving the M(3,p) Virasoro minimal models that has been considered 
previously, using results from the literature. We exhibit the Hamiltonians for the first few cases of 
these. 



I. INTRODUCTION 



Overview 



In recent years there has been a considerable growth 
of interest in topological phases of matter (see e.g. Ref. 
0, [|[ for reviews). A topological phase of matter can 
be defined as an equilibrium zero-temperature phase of 
a quantum system that has a gap in its energy spectrum 
above the ground state, as far as local excitations far 
from any boundary are concerned. Topological phases 
may have non-trivial properties (also called topological) 
that are unchanged by any (local) perturbation of the 
(local) Hamiltonian unless it causes the system to pass 
through a phase transition. Model examples of systems 
(consisting of a space of states, and a Hamiltonian) ly- 
ing in a particular topological phase provide existence 
proofs for the phase. Due to the gap in the spectrum, 
any sufficiently small perturbation of the Hamiltonian 
by a local term will not destroy the gap. Hence, the 
solutions (even partial solutions) of such model exam- 
ples provides paradigms for those phases. Among the 
earliest examples of topological phases and of the use of 
partially soluble models is the fractional quantum Hall 
effect (FQHE) 0,13. 

Historically, a standard approach has been the use of 
trial (or variational) states for the ground and some of the 
excited states, even in the absence of a soluble Hamilto- 
nian. In the FQHE, examples of such states are usually 
described by a wavefunction that may readily be writ- 
ten down for arbitrary size (particle number) in terms of 
multiplying some factors, performing integrals or sums, 
or other operations. It is then hoped that these trial or 
"nice" wavefunctions have the topological properties of 



a topological phase. In some cases, a Hamiltonian (or 
rather, a family of Hamiltonians, one for each system 
size) is known for which the ground state is also known 
exactly for each size of the system. As generic states are 
hard to specify concisely, such a state must have some 
special structure that enables it to be specified in a sim- 
ple way; thus it too must be a nice function. 

Various approaches to obtaining trial wavefunctions 
have been used in the past. In the FQHE, after Laugh- 
lin's wavefunctions Q, there were generalizations using 
the hierarchy and composite particle approaches 043) 
though some of these schemes possess considerable lee- 
way when it comes to specifying definite functions. In 
this paper we will be concerned with another approach 
that developed from Laughlin's states. It stemmed from 
the observation that the algebraic behavior of the func- 
tions, which vanish as a power of separation as any two 
particles approach the same point, implies that there is a 
"pseudopotential" interaction Hamiltonian acting in the 
space of lowest-Landau level (LLL) states that annihi- 
lates both the Laughlin ground state and also the states 
containing any number of fractionally-charged quasihole 
excitations [4|. Thus the "nice" wavefunctions are actu- 
ally zero-energy eigenstates of a "special" Hamiltonian, 
and also exhaust all such states. The special Hamilto- 
nian consists of projection operators acting within the 
LLL; the coefficients of the projectors are positive, so 
the Hamiltonian is positive (semidefinite) as an opera- 
tor within the LLL. Moreover, the special Hamiltonian 
is usually translationally and rotationally invariant. It 
is also local (short range), because it involves the low- 
degree behavior of the wavefunction in the differences of 
particle coordinates as they tend to zero. 

Another approach to obtain trial wavefunctions was 
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developed by Moore and Read (MR) [7], who ob- 
served that the Laughlin wavefunctions, and some other 
previously-known trial wavefunctions, are chiral correla- 
tion functions (so-called conformal blocks) in some two- 
dimensional conformal field theory (CFT) [1,0. They 
extended this observation to propose some further trial 
states, and the method used amounts to a technique for 
the construction of an unlimited number of trial ground 
and associated quasihole wavefunctions. Subsequently, 
it was observed fiol . that the main example in MR, 
the MR or Pfaffian state, is itself the zero-energy ground 
state of a special (i.e. projection operator) Hamiltonian. 
This Hamiltonian is similar to those for the Laughlin 
state, in that it involves the way in which the function 
vanishes as any three particle coordinates come to the 
same point, and so consists of a three-body interaction. 
Using this Hamiltonian, explicit wavefunctions for any 
number of quasiholes were constructed [TH [ill- Other ex- 
amples with similar structure have been found, including 
notably a sequence of further distinct topological phases 
in which the special Hamiltonians contain k + 1-body 
interactions, where k = 3, 4, 5, ...[13]; the preceding 
examples fit into the sequence as the cases k = 1, 2. 

Of the sets of ground and quasihole trial LLL wave- 
functions that have been proposed, many of those that 
can be expressed as functions, without use of integration 
over dummy variables (such as use of projection opera- 
tors, possibly from a larger Hilbert space such as those 
including higher Landau levels), can be viewed as con- 
formal blocks in a (known) CFT. Further, to our knowl- 
edge, all such sets of functions, that appear acceptable 
as a space of ground and quasihole states in a topological 
phase and are precisely the zero-energy states of a LLL 
special Hamiltonian, are conformal blocks. It should be 
emphasized that even when all the zero-energy states of 
a special Hamiltonian can be found, the non-zero energy 
eigenvalues generally cannot, and one does not know di- 
rectly whether there is a gap in the spectrum of bulk exci- 
tations over the ground state in the infinite size limit. In- 
deed, it has been argued that when the wavefunctions are 
conformal blocks of either a non-rational or non-unitary 
CFT, any such Hamiltonian must be gapless in the limit 
[Til [l5j . Consequently, in these cases the wavefunctions 
do not represent the behavior of a topological phase, but 
may be at a critical point. 

In this paper, the main goal is to provide tools to study 
systematically the short-distance behavior of wavefunc- 
tions (and related entanglement properties), use this to 
find special Hamiltonians, and analyze the spaces of all 
zero-energy states of those Hamiltonians. The proper- 
ties of the wavefunctions in these spaces lends itself to 
an algebraic interpretation in general, and leads almost 
inevitably to the use of MR constructions from a CFT to 
construct and analyze examples. In giving an overview 
of these ideas and results to be presented in this paper, 
we will now divide them under the following headings. 



1. Zero-energy subspaces for a given Hamiltonian 

LLL wavefunctions for bosons are essentially symmet- 
ric polynomials in the complex coordinates Zj = Xj + iy.j 
for the jth particle with Cartesian coordinates Xj , yj in 
the plane. From here on, when we refer to functions, it is 
such symmetric polynomials that we mean. The phrase 
"short-distance behavior" means a Taylor expansion in 
some selected particle coordinates as they all tend to 
some value; hence we usually organize the spaces by de- 
gree of the polynomials, beginning with the lowest degree. 
(We concentrate on bosons; wavefunctions for fermions 
may be obtained by multiplying by Jli<j( z i — z j)t an< ^ 
there are related Hamiltonians.) Given a translationally- 
and rotationally-invariant projection operator Hamilto- 
nian that consists of interactions among at most N max 
particles, one can in principle determine all of the zero- 
energy states for each number of particles TV — 1, 2, ... in 
the plane. We call the spaces of these states the zero- 
energy (sub-) spaces of states or wavefunctions; they are 
subspaces of the space of all LLL states for given N, and 
usually only of interest when they are proper subspaces. 
We also describe the wavefunctions as having "allowed" 
behavior (one might also describe them as "satisfying" 
the constraints that they are annihilated by the projec- 
tors in the special Hamiltonian). In the literature, these 
are often called spaces of edge excitations, though they 
can only be considered as true edge excitations in the 
case when there is a gap for bulk excitations [H, Il7j . 
(One may also consider other geometries, including the 
sphere with a given number of flux quanta emerging from 
it, which lead to related results. For a gapped bulk phase 
on the sphere, the zero-energy states may be viewed as 
quasihole states.) In the body of this paper, we carry out 
the analysis of the zero-energy subspaces for the plane 
and the sphere explicitly in some examples. 

2. Entanglement subspaces for a given ground state 

In analyzing the zero-energy states, the spaces for 
N < N max play a special role. If these are known, then 
the zero-energy wavefunctions for N > iV max particles 
must be found by some process of "patching together" the 
zero-energy wavefunctions for sets of Na < -/V max parti- 
cles. Put another way, if the former are decomposed into 
wavefunctions for a subset of Na < N max particles and 
for the remainder, the wavefunctions for the Na particles 
must lie in the zero-energy subspace. Indeed, the same 
is true for the decomposition for any value of Na < N. 
For the "ground state" in the zero-energy subspace for 
any N, that is the (usually unique) wavefunction of low- 
est total degree in the z^s (lowest angular momentum), 
we term the space of functions that actually occur when 
it is decomposed into functions of Na, times functions 
of the remaining N — Na , coordinates an entanglement 
subspace, because of its relation (discussed below) to en- 
tanglement properties of a ground state wavefunction. 
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Prom the preceding remarks, it is clear that the entangle- 
ment subspaces for any Na (for the ground state at any 
N > Na) must lie within the zero-energy subspace for 
N = Na particles. [One could of course do the same for 
any zero-energy state; we sometimes refer to these also as 
entanglement subspaces.] In the body of this paper, we 
analyze the relation of the zero-energy and entanglement 
subspaces, in some examples. 



3. Inverse problem: finding a special Hamiltonian for given 
zero-energy subspaces 

If there are some polynomials "missing" from the en- 
tanglement subspace for some value of Na, and some 
degree d, for all values of the total number N, then a 
projection operator onto that function can be added to 
the Hamiltonian, and it will still annihilate the wavefunc- 
tions. More precisely, the projection must be onto some 
part of the orthogonal complement of the "allowed" func- 
tions. (The orthogonal complement and the construction 
of projection operators involve the use of the quantum- 
mechanical inner product on the wavefunctions of our 
system.) The special Hamiltonian with which we began 
must itself consist of such projectors, and one would like 
to use one with a minimal number of terms. Then we also 
have the inverse problem: if we are given some spaces of 
wavefunctions, for all particle numbers N, and we wish 
to obtain a special Hamiltonian for which these are pre- 
cisely the zero-energy subspaces, then such an analysis 
of the entanglement subspaces will ultimately lead to the 
solution, if one exists. A necessary condition for existence 
is that the spaces of functions should obey the compat- 
ibility condition mentioned in the preceding paragraph. 
That is, all the entanglement subspaces for Na particles 
in any iV-particle function in the given space (Na < N) 
should be a subspace of the given space of functions for 
Na particles. 



4- Spaces of amplitudes from a CFT 

In particular, as spaces of functions, we can begin with 
conformal blocks in a MR construction. This produces 
spaces of functions, which we may also term "spaces of 
amplitudes" , that play the role of the zero-energy sub- 
spaces; the spaces depend on the CFT used. In these 
constructions, the necessary compatibility conditions are 
automatically satisfied, because the short-distance anal- 
ysis of entanglement subspaces is related to the operator- 
product expansion in CFT, which has such behavior f]~3j j . 
We may attempt to find a corresponding special Hamil- 
tonian for which such spaces of amplitudes are the zero- 
energy subspaces. 



5. Entanglement 

The analysis so far has led to issues of entanglement. 
For entanglement in general in a many-particle state, one 
effectively divides the system into one part with Na par- 
ticles, and the remainder with Nb = N — Na particles. 
One may trace over the states of the Nb particles to 
obtain a reduced density matrix, and then analyze its 
eigenvectors and eigenvalues QJ, HU. Thus essentially, 
the space in which the eigenvectors lie (the entanglement 
subspace) is that spanned by the wavefunctions of the 
Na particles that occur in the ground state, just as in 
the preceding discussion (we will discuss different ways 
of defining the bipartition of the system in more depth 
later). If the partition respects rotation invariance about 
the origin, then the states can be classified according to 
their angular momentum, which is the same as their de- 
gree as polynomials. In some ways of defining the bi- 
partition to obtain entanglement, a range of values of 
Na must be considered also. Several ways of biparti- 
tioning for LLL states have been considered in the litera- 
ture; for our purposes, which relate to special Hamiltoni- 
ans, the most appropriate is particle partition [20], which 
preserves translational, as well as rotational, symmetry. 
(Many studies of entanglement use the sphere; in that ge- 
ometry particle partition has the full rotation symmetry.) 
When studying entanglement, one is often interested in 
universal behavior as Na (or its "typical" value, if there 
is a range) tends to infinity (either subsequent to, or at 
the same time as, N — > oo). 

An important idea is that in LLL trial states for a topo- 
logical phase in the FQHE (in fact, for those obtained as 
conformal blocks/amplitudes from a MR construction), 
the multiplicities of the states in the decomposition at 
sufficiently large angular momentum and Na — that is, 
the dimensions of the entanglement subspaces — approach 
those of the edge states [21[, which can sometimes be 
identified as the zero energy states of a special Hamilto- 
nian. For states given by the MR construction as confor- 
mal blocks, these multiplicities (in the limit) are those of 
the vacuum sector of the CFT. More generally, in that 
construction, the entanglement subspaces must lie within 
the spaces of amplitudes, because of the compatibility 
conditions already mentioned. We will show, by general 
algebraic arguments, that the dimension of the space of 
amplitudes obtained by the MR construction does agree 
with the vacuum sector of the CFT used as N — > oo (in 
an appropriate way). We also show that, in the limit 
Nb — > oo with Na fixed, the entanglement subspaces 
contain all the same functions as those in the spaces of 
amplitudes for N — Na, at any finite value of Na as 
well, provided the correct CFT is used to construct am- 
plitudes. Hence in particular, their dimensions are equal. 
This is a finite size (Na finite, though Nb went to infin- 
ity) version of a bulk-edge correspondence, if the spaces 
of amplitudes are viewed as edge states. The "correct" 
CFT is one without any "singular" vectors, or alterna- 
tively can be obtained from the given CFT by a process 
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of setting the singular vectors to zero, to form the quo- 
tient space [H (Singular vectors can only occur in 
non-unitary CFTs.) The equality of the entanglement 
subspaces with the spaces of amplitudes does not hold 
for general Na as Nb — > oo if there are singular vectors 
in the CFT that was used to find the amplitudes, even 
though the ground state wavefunctions, that result from 
the construction using the CFT instead of its quotient, 
are the same. We give examples of this phenomenon. 



6. Algebraic techniques 

We develop an algebraic language to aid in all of this 
analysis. The space of all symmetric polynomials in Na 
variables is a vector space, and can be viewed as the 
Hilbert space of LLL states of bosons. But it can also be 
viewed as an algebra, because not only are sums and com- 
plex multiples of symmetric polynomials again symmetric 
polynomials in the same variables, but so are products of 
them. Further, this algebra acts in an obvious way, by 
multiplication, on the vector space of symmetric poly- 
nomials, which can then also be viewed as a "module" 
over (or in less fancy language, a representation of) this 
algebra. 

It is very convenient if the zero-energy subspaces, 
which are subspaces of symmetric polynomials, are them- 
selves modules over the symmetric polynomials — that is, 
if the subspace is closed under multiplication by sym- 
metric polynomials. That is because, as we will see when 
explicitly finding allowed functions, this property allows 
us to combine factors for different subsets of the parti- 
cles that satisfy the conditions on allowed functions sep- 
arately, secure in the knowledge that they will also do so 
when combined as N particle wavefunctions. This prop- 
erty holds in previously known examples, and in those 
we analyze here. If the special Hamiltonian consisted of 
an arbitrary set of projectors (i.e. if the allowed functions 
did not form a module over the symmetric polynomials), 
it would be much harder to construct the spaces of zero- 
energy many-particle states, and it would not be clear if 
any such states even exist. We further require that our 
Hamiltonians be short-range interactions, and that they 
be translationally and rotationally invariant. Rotational 
invariance implies that angular momentum, or the degree 
of the polynomials, be a good quantum number; we may 
classify functions by degree. Translational invariance im- 
plies that the allowed behavior is mapped into itself by 
the application of the center of mass variable cx JV z% 
and J2i d/dzi (acting on the polynomials). 

Consequently, when searching for special Hamiltoni- 
ans, we now have an algebraic problem, involving the 
study of translation-invariant modules for the algebra 
of symmetric polynomials. While this problem can be 
solved completely for the case of Na — 2 particles (see 
Section II CI below), it is much more complex for more 
than two particles. To make progress, motivated by the 
MR construction using CFT, we can also enlarge the al- 



gebra of symmetric polynomials by introducing further 
operators which form part of the action of the Virasoro 
algebra on the wavefunctions. This algebraic problem 
can also be studied, without further use of CFT. Be- 
cause the requirement of being a module over the larger 
algebra is more stringent, there are fewer solutions to be 
found, and they potentially arise from conformal blocks. 
Conversely, given such a Hamiltonian, we may find its 
zero-energy states, as we will do in several examples. 

Further motivation for confining the search to spaces 
with such an (extended) module structure can be ob- 
tained from arguments in Ref. [l5j]. It was argued there 
that for a special Hamiltonian that has a gap in the bulk 
spectrum (and no other assumptions) , the space of zero 
energy (or edge) states in the limit of large sizes forms 
a chiral conformal field theory, with an action of both 
the Virasoro and U(l) current algebras. However, it is 
not clear if all the operators that generate these algebras 
(or of the subalgebras that we use in this paper) must 
have the one-body form that we encounter below, even 
though certain of them, which relate to translations and 
rotations of all the particles, do. (The arguments of Ref. 

also lead us to expect that the action of the current 
and Virasoro generators must be local in position space; 
it may be that this further condition implies that the one- 
body form must hold.) Nonetheless, such a form would 
seem natural, and the assumption that there is this mod- 
ule structure, which we will make even for cases that we 
don't expect to be gapped in the bulk, and in finite size, 
will ensure that when the bulk is gapped (the cases of 
most interest to us anyway), this physics will emerge. 

In the spaces of amplitudes produced by a MR 
construction, the entanglement subspaces automati- 
cally form a module over the symmetric polynomials — 
corresponding to the presence of a U(l) current algebra 
in the CFT — and over the part of the Virasoro algebra 
mentioned above, and perhaps also over some larger al- 
gebra, depending on the CFT considered. In this set-up, 
the "missing" polynomials (corresponding to the orthog- 
onal complement of an entanglement subspace) are re- 
lated to a finite-size version of the presence of singular 
(also called null) vectors in the underlying Verma mod- 
ules of the CFT, a topic familiar in CFT. If we under- 
stood these orthogonal complement spaces, we would be 
able to find projection operators that could be included 
in a special Hamiltonian whose zero-energy spaces con- 
tain the given spaces of amplitudes. We would then be 
left with the problem of finding a minimal Hamiltonian 
such that the latter two sets of spaces actually coincide. 



7. Superconformal and M(3,p) examples 

To illustrate the preceding ideas in a concrete manner, 
we introduce a family of examples. The original motiva- 
tion for this work was to find examples of special Hamil- 
tonians for which the zero-energy states are precisely the 
spaces of amplitudes obtained from a unitary rational 
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CFT; a guiding example was the tricritical Ising state 
14 1 for which the ground state wavefunction is known 
22]. Unfortunately we have not succeeded in finding a 



special Hamiltonian for this example. But Ref. [22| gives 
a family of ground state wavefunctions, one for each of a 
continuous family of special Hamiltonians, and those are 
labeled by one parameter c and involve three-body inter- 
actions only. These wavefunctions are obtained from the 
MR construction using the N = 1 superconformal field 
theories, and include the tricritical Ising ground state at 
one special value of c, the central charge. (In general, 
these SCFTs are neither unitary nor rational.) The fill- 
ing factor in these states is v = 1/3. 

In Section |TI] of this paper we analyze the zero-energy 
spaces for this family of special Hamiltonians (with no 
use of CFT), and find their dimensions. We also show 
that this family of models based on Af — 1 SCFT, which 
are expected to have gapless bulk excited states in the 
limit, are adjacent to a MR phase of bosons at v = 1/3, 
consistent with their being on a phase boundary. In Sec- 
tion ing we show that (except for c = 0, which has larger 
degeneracy) in the case of the plane geometry these di- 
mensions agree in the infinite size limit with those of the 
generic Af = 1 superconformal field theory (SCFT) con- 
struction. The amplitudes constructed from the SCFT 
are guaranteed to be zero-energy states, and in fact we 
can eventually prove that (for c ^ 0) the zero-energy sub- 
spaces for any finite N are identical with the spaces of 
amplitudes of that SCFT for the same N (that is, not 
only the ground state wavefunction, but all wavefunc- 
tions in these spaces). Hence these spaces are also equal 
(at generic values of c, as we discuss in a moment) to the 
entanglement subspaces at Na equal to this value N (as 
Nb —> oo); the functions in, and the dimensions of, these 
spaces are then all known. These arguments are very 
general and can be applied to many other MR construc- 
tions, whenever the zero-energy states can be constructed 
and the dimensions of the spaces agree with those of the 
vacuum of the CFT in the large-size limit. 

At the special point corresponding to the tricritical 
Ising model, these generic dimensions are strictly larger 
than those for the spaces of amplitudes in the latter. This 
is because at this value of c, the modules of the generic 
Af = 1 SCFT (which will be precisely defined later) con- 
tain singular vectors. The quotient space, obtained by 
setting the singular vectors to zero (or "modding them 
out"), yields a module that has no remaining singular 
vectors. The same phenomenon occurs at a countable 
infinity of other "rational" values of c, and is connected 
with the Af — 1 SCFT minimal models. In general, the 
SCFT minimal models (like the Virasoro minimal mod- 
els) are rational but not unitary, however, there is a sub- 
set that are also unitary; among the latter, the tricritical 
Ising theory is the non-trivial theory of lowest c. 

These facts have several consequences. First, at these 
rational values of c, the entanglement subspaces of the 
ground state, which as we will show always agree with 
the spaces of amplitudes in the MR construction using 



the quotient or minimal model SCFT, have lower dimen- 
sions than those at generic c. This shows, as mentioned 
above, that the correspondence of dimensions of the en- 
tanglement and zero-energy spaces does not always hold. 
Hence, as c varies, some pseudoenergies in the entangle- 
ment spectrum must move off to infinity at these values of 
c, and this will occur even in finite size (sufficiently large, 
though how large will depend on the c value involved). 

Second, at the rational values of c, the "missing" poly- 
nomials in the entanglement subspaces for some values 
of Na allow us to add corresponding terms to the spe- 
cial Hamiltonian which then still annihilates the ground 
state. These terms depend on c and do not annihilate 
the ground state at generic c. 

In order to determine, for a particular rational value 
of c, whether addition of such terms in the Hamiltonian 
produces zero-energy states that coincide with the spaces 
of amplitudes in the MR construction using the quotient 
SCFT at the same c, one must analyze the zero-energy 
states, and then if they are not the same as the spaces of 
amplitudes, add further terms to the Hamiltonian. This 
seems difficult to carry through in general. However, a 
remarkable paper by Feigin, Jimbo, and Miwa (FJM) 
[23^ (whose algebraic approach is very similar to ours, 
though without the physical applications in view here), 
shows in effect that for one rational value of c, there 
is a three-body special Hamiltonian whose zero-energy 
subspaces are precisely the spaces of amplitudes, for any 
particle number TV. In fact, their analysis applies to a 
sequence of MR constructions using Virasoro minimal 
models M(3,p) [p is a positive integer, and not divisi- 
ble by 3; p = 4 gives the MR state, and M(3, 5) is the 
Gaffnian [2J]]. These spaces of amplitudes have been 
discussed previously in relation to the FQHE, but the 
existence of a special Hamiltonian in these models was 
apparently overlooked (except for p = 4, 5, for which 
the Hamiltonians were known). We give the Hamiltoni- 
ans for several of these in explicit form, and find all their 
zero-energy wavefunctions for all p for both the plane and 
sphere geometries; the multiplicities of these reproduce 
results in FJM. 



8. Structure of the paper 

The structure of this paper is as follows. In the re- 
mainder of this section, we review some of the earlier 
work that is relevant for our paper, and describe the is- 
sues in greater detail. The basic issues are dual: spe- 
cial Hamiltonians and subspaces of "nice" trial functions. 
This duality leads to a "chicken and egg" problem, in 
that a special Hamiltonian (the chicken) produces spaces 
of nice zero-energy states (the eggs) , while in turn spaces 
of states that satisfy the compatibility conditions can be 
used to produce Hamiltonians that annihilate them. It is 
not obvious with which of these two types of objects one 
should begin. But as both types of analysis hinge on the 
entanglement subspaces, which correspond to (possible) 
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terms in a special Hamiltonian, we begin with those. 

In Section [Til we consider the wavefunctions deter- 
mined by the (continuous family of) "vanishing condi- 
tions" on the relative variables for any three particles, 
which is related to the TV = 1 SCFT. The methods 
make no use of CFT, and give the complete, linearly- 
independent states for both the plane and sphere geome- 
tries. The count of states on the plane is interpreted later. 
We also argue that this model sits at a phase boundary 
adjacent to a MR phase. 

In Section Mil we begin to use methods from CFT. 
We describe features of the MR construction in the op- 
erator language of vertex operator algebras (VOAs). We 
prove a Theorem that implies that the functions obtained 
from the MR construction span a space isomorphic to the 
VOA. We show that the entanglement subspaces of the 
ground state in particle partition span all the amplitudes 
in the limit as Nb — >• oo if the CFT is non-degenerate, 
but in degenerate cases span a smaller space. We show 
that in CFT-related examples in which the count of zero- 
energy states of a special Hamiltonian can be carried out, 
these numbers and the wavefunctions themselves agree 
with the amplitudes even in finite size TV. 

Finally, in Section IIV1 we show in an example how 
when the entanglement subspaces are smaller than the 
space of amplitudes of a CFT (or zero-energy states of a 
special Hamiltonian), because the CFT contains singu- 
lar vectors, terms can sometimes be added to the spe- 
cial Hamiltonian such that agreement is restored with 
the quotient CFT. The examples in this section use the 
M(3,p) minimal models. 

In Appendix [Al the proof of Theorem 1 is given. 



B. Decomposition of trial wavefunctions 

We consider the example of the Laughlin-Jastrow 
wavefunction Q 

*T J (z U ...,Z N )=]J(z i -Z j ) Q , (1.1) 

i<j 

where we have dropped the ubiquitous factor that is a 
function of \zi\ 2 (which is Gaussian for the plane [3j, ra- 
tional for the sphere [ll|) as they do not enter the argu- 
ment here. The degree in each Zi is (the number of 
flux quanta covered by the particles), 

N 4 , = QN-Q. (1.2) 

If we include quasihole factors Oil^j — w k) f° r each of 
a set of quasiholes at positions Wk (k = 1, . . . , n), then 
we know that each quasihole factor can be expanded in 
elementary symmetric polynomials 

N 

Y[(zi -w k ) = J2 e m (-w k ) N - m , (1.3) 

i m—0 



where 

f-m ^ Z{ 1 Zi 2 • • • Zi m (1.4) 

il<i2<---i m 

(to = 1, . . . , TV) are the elementary symmetric polyno- 
mials in TV variables, and for convenience we also define 
e o by eo = 1. The elementary symmetric polynomials 
generate the algebra of all symmetric polynomials in TV, 
and do so freely (there are no linear relations between 
products of elementary symmetric polynomials). Then 
expanding the product of quasihole factors in powers of 
the Wfc's produces symmetric polynomials, and by using 
sufficiently many quasiholes we can eventually produce 
all the symmetric polynomials. [As the functions are 
symmetric under permutations of the w k s, this expan- 
sion can be viewed as a double expansion in symmetric 
polynomials in the two sets of variables]. These func- 
tions which contain the Laughlin-Jastrow factor times 
arbitrary symmetric polynomials can be viewed as rep- 
resenting the edge excitations of the Laughlin disk. 

The Laughlin wavefunction itself has a somewhat sim- 
ilar form as the quasihole factors. Thus if we di- 
vide the particles into two groups, say {1, . . . , N A } and 
{N A + 1, . . . , TV} (and let N B = TV - N A ), then we may 
write the function in the form 

*L(Z1) • • • ,*n) = y]lf)Ar(zi,.^,Zff A ) 
r 

X1p Br (zN A +l,---,Z N ), (1.5) 

where r is an arbitrary index and the functions ipAr and 
ipBr are symmetric in their arguments. If ^ l were re- 
placed by a generic state (symmetric polynomial), the 
functions ipAr obtained in this way would span the set 
of symmetric polynomials times V'l in Na variables, at 
least in low degree (clearly the degree of ip A r must be 
less than TVTV^/2, the degree of the state considered). It 
is natural to wonder if this is the case for the Laugh- 
lin state. There are several variants of this question (in 
each case, initially with Na fixed): (i) in one version, 
we consider all polynomials obtained as Nb increases to 
infinity; (ii) in a second, Nb is instead fixed, and it fol- 
lows that the highest degree in any z% with i G A is 
bounded by [version (i) corresponds to the limit of 
this as Nb — > oo]. In one further version, (iii) we require 
that Nb be fixed, and also that in the functions ipAr the 
degree in each in {1,2,..., Na} be less than or equal 
to some bound less than TV^, while the degree of ip B m 
in each Zi in {Na + 1, . . . , N} must be strictly greater 
than the same bound, similarly. In this case, an explicit 
symmetrization over permutations of the coordinates is 
required in order to yield a valid decomposition of 
in eq. (jl.5l) (it is not required in the others preceding). 
In each version, one may finally consider different Na, 
and take the limits as Na — > oo. These decompositions 
can be applied in the same way to other trial wavefunc- 
tions also, and the limits can be studied as there is such 
a wavefunction for each TV (or for each TV divisible by 
some k > 0). 
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These decompositions of the Laughlin wavefunction 
are equivalent to definitions of entanglement that have 
been considered recently. For entanglement, we assume 
that the Hilbert space % has a tensor product form, 

U = Ua®U b (1.6) 

at least, within each sector of fixed Na- If we represent 
vectors in each space by wave functions, then we can de- 
compose a function ip accordingly as 

1p = ^ ^At^Bt-, (1-7) 
r 

where ipAm (ipBm) belongs to Ha ("Hb)- The question 
is to characterize the degree and nature of the entan- 
glement, that is the extent to which ip is not a tensor 
product. (This is rather too naive for systems of many 
identical particles, which because of symmetrization do 
not have a product form. However, once some particles 
have been assigned as part A and others as B, the overall 
symmetry of the state is broken, and the product form 
applies.) One version for QH states corresponds to fur- 
ther restricting the single-particle LLL Hilbert space ac- 
cording to the angular momentum of each particle, which 
corresponds to the degree in each Z{, so that some values 
appear only in part A, and the remainder in part B. This 
"orbital" entanglement [2(| [2l[ corresponds to the third 
version above. The entanglement subspaces spanned by 
the tpAr for some trial wavefunctions were found to be 
smaller than the spaces of all symmetric polynomials of 
the correct degree — this is essentially the phenomenon 
of the "entanglement gap" discovered in Ref. [2l]L_ Other 
studies have considered "particle" entanglement [2fJ , cor- 
responding to the second version above. (A more formal 
extended definition of particle partition, in particular, 
was given in Ref. [25].) It has been conjectured that 
the dimensions of the entanglement subspaces in parti- 
cle partition on the sphere agree with the counting of 
quasihole states for N = Na, N^, provided N B is suf- 
ficiently large [2^, however this does not seem to have 
been proved to our knowledge. Another version of the 
entanglement spectrum is obtained by considering a bi- 
partition in position space — the real-space entanglement 
spectrum [25|, [27h[29| , which we will not consider further 
here. The limit that has usually been considered for or- 
bital and for real-space entanglement is that in which 
N — s- oo with the ratio Na/N fixed. The first version 
above (Nb — > oo with Na fixed) has not been considered 
as much. 

As mentioned above, in Section lHTl of this paper, we ob- 
tain detailed statements about the entanglement spaces 
in particle (and also real-space) entanglement, in the first 
version above, the limit N — > oo with Na fixed. We 
believe this is the simplest, or most basic, limit of all. 
That is, we can show that after taking Nb — > oo, the 
functions ipAr span a space isomorphic to the space of 
amplitudes in Na variables, where the CFT used to con- 
struct the amplitudes is a quotient of the CFT used in 



the construction, obtained by setting any singular vectors 
to zero. Moreover, as Na — > oo, the space of amplitudes 
in a MR construction becomes isomorphic to the space 
of states in the vacuum representation of the same CFT. 
We note that there has been recent progress in establish- 
ing such relations for real-space and particle partitions 
of trial wavefunctions [29j . with which the present re- 
sult has some overlap. However, the arguments presented 
there made use of generalized screening behavior to ob- 
tain results about the inner products of the wavefunc- 
tions, which yields much more information about the en- 
tanglement spectrum (not only dimensions of subspaces) 
than we obtain here. Such behavior is not expected to 
hold for trial wavefunctions obtained from non-unitary 
or non-rational CFTs. By contrast, the method used in 
the present paper is purely algebraic in that the quan- 
tum mechanical inner product on the wavefunctions is 
not used, and while it gives less information, it works for 
all cases, including non-unitary and non-rational ones. 

C. Subspaces and Hamiltonians 

As we know, the entanglement subspaces are the key 
to searching for a special Hamiltonian for a given set 
of wavefunctions. Several works have addressed this 
problem, beginning from Ref. [4[. Recently, a start has 
been made on systematically studying the possible spe- 
cial Hamiltonians for more than two-body interactions, 
and finding their zero-energy many-particle states [30| . 

It will be useful to describe the spaces of symmetric 
polynomials in N variables in more detail. Because the 
algebra of symmetric polynomials is freely generated by 
the elementary symmetric polynomials (or, in fact by 
some other alternative sets of N generators), the num- 
ber of symmetric polynomials in each degree can be found 
easily. Generally, for a vector space V of symmetric poly- 
nomials that have multiplicities (dimensions) bd in each 
degree d, we can define the character, 

oo 

ch g A = Y. bd( l d ( L8 ) 

d=0 

where q is an indeterminate (the series is a formal power 
series, which means that its convergence is not of in- 
terest). For the space An consisting of all symmetric 
polynomials in N variables, the character is 

N 1 1 

where we have introduced the compact and standard no- 
tation, 

m 

(q) m = Y[( 1 -1 n y ( L1 °) 

n=l 

Because the Hamiltonian should be translationally (and 
rotationally) invariant, it is useful to decompose the pos- 
sible behavior for N particles into the center of mass and 
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some "relative" or "internal" wavef unction. The center 
of mass coordinate is Z = ^2iZi/N — e±/N , and the cen- 
ter of mass wavefunctions are spanned by Z m , m = 0, 
1, 2, .... The center of mass angular momentum can 
be raised by multiplication by Z, and lowered by apply- 
ing d/dZ = J2i^/dzi. (In these expressions, the sums 
over particles are from i = 1 to N. We can treat the 
derivatives as acting on the polynomials part of the wave- 
function only [3lj). For the internal part, coordinates 
can be defined, and a space of functions of these coor- 
dinates orthogonal to the center of mass variable can be 
constructed. However, there is no particularly natural 
definition for these. In any case, the functions in the 
internal variables must be annihilated by d/dZ, which 
means that the polynomials are translation invariant — 
they are functions of differences Zi — Zj only. The space of 
general wavefunctions is thus spanned by Z m times poly- 
nomials of the differences of coordinates. (The space of 
translationally-invariant symmetric polynomials can be 
freely generated by a modified set of N — 1 elementary 
symmetric polynomials, given in Ref. 30].) For two par- 
ticles, the latter have the simple form 



(i.ii) 



where again m' = 0, 1, . . . , and now m' must be even for 
bosons (symmetric polynomials) and odd for fermions. 
There is only one polynomial for each even (resp., odd) 
degree. This is not typical; for N > 2, the number 
of linearly-independent polynomials in each degree in- 
creases gradually. The space of polynomials in Z has 
character 1/(1 — q), and so the character of the space of 
wavefunctions for the internal motion is 



n 



i 



(I-?™)' 



(1.12) 



For example, for N = 3 particles, we have 



1 



i+q 2 +q 3 +q' l +q 5 +2q 6 +q 7 +2q 8 - 



(l-g2)(l-g3) 

(1.13) 

and the coefficients have the pseudoperiodic property 
that if the degree increases by 6, the coefficient increases 
by 1. 

Following earlier examples, we may attempt to de- 
scribe some spaces of polynomials by specifying the be- 
havior that is "allowed" for members of the space as some 
numbers Na of the particles come together; these are 
sometimes termed "vanishing conditions" . As discussed 
above, we choose to require the spaces to be modules over 
the algebra of symmetric polynomials. It is not clear if 
this property is really necessary. Physically, multiplying 
by a symmetric polynomial generates an excitation of the 
charge sector only. 

For Na — 2, the translation- invariant modules over the 
symmetric polynomials can be completely determined. 
There is only one generator, (z\ — Z2) 2 , so any such mod- 
ule is spanned by the even powers {z\ — Z2)™ 1 as above, 



for m' =m, m + 2, . . . , times a power Z M , M — 0, 1, 
.... Thus the distinct translation-invariant modules are 
labeled by m — 0, 2, . . . , and the elements of each such 
module can be generated by (z% — z^)™ 1 by multiplication 
by symmetric polynomials. These correspond to the dis- 
tinct special or pseudopotential Hamiltonians for bosons 
that produce the Laughlin ground, quasihole and edge 
states as zero-energy states [J] (it is similar for fermions, 
if we consider modules consisting of antisymmetric poly- 
nomials). Cases that include conditions on more than 
two particles are much richer. They include the MR 
state for bosons at v = 1, which arises from requiring the 
polynomials to vanish if any three coordinates coincide. 
This corresponds to the degree term in the above series, 
which represents the dimension of the relative-coordinate 
space of symmetric polynomials. Because there is no rel- 
ative term at degree 1, the functions vanish quadratically 
as the three coordinates come together (in any fashion, 
e.g. first two, and then the third, or with any fixed ra- 
tio of pairwise separations) . The more-recently proposed 
Gaffnian ground state, and its edge and quasihole states, 
arises if we require the functions to vanish as degree 3 
when three coordinates approach the same point [24j . 
Both of these spaces of allowed behavior form a module 
over the symmetric polynomials. 

In the examples considered in Ref. [3(J, the allowed 
behavior was defined by requiring that all relative mo- 
tion with degree less than some bound was forbidden, 
and all greater was allowed. We note that this automati- 
cally produces a module over the symmetric polynomials. 
However, not all those examples led to truly "nice" be- 
havior, with a unique ground state (lowest total degree 
polynomials) for N particles, and moderate growth of 
the number of allowed or zero-energy states at higher de- 
gree. Based empirically on the evidence of past cases, 
such behavior would be highly desirable when searching 
for new model QH phases of matter. Other than the re- 



cent Gaffnian [24| , no such nice spaces of functions were 



identified in that work that had not been found previ- 
ously. In further work (32|, some four-body interactions 
were studied, and related to a MR construction with S3 
CFTs. The methods developed in the present paper can 
be applied to that case also. 

At the higher degrees, for example degree 6 or larger 
than 7 for the three-particle case above, there are more 
possibilities. We may choose to allow, for three particles, 
any degree greater than 6, and only a one-dimensional 
subspace of the space in degree 6. The choice of a sub- 
space means that there is now a family of spaces of func- 
tions (a one-parameter family in this example), or a fam- 
ily of Hamiltonians, each producing one of the spaces. 
We will study this example in detail in this paper. We 
emphasize that the possibility of such a parameter or pa- 
rameters was not included in the original versions of the 
"thin-torus-limit" . "pattern-of-zeroes" , and "root state" 
approaches [33l - l39j . But the objects under study are vec- 
tor spaces, which allow use of such linear combinations 
(when the dimension is greater than one), not discrete 
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objects such as sequences of integers. The latter cannot 
represent general vectors in a vector space, though they 
can represent basis vectors, such as the monomial ba- 
sis for many-particle states in the lowest LL. We should 
add that later versions of some of these approaches have 
gradually included more terms in the Hamiltonian, corre- 
sponding to subleading terms in the thin-torus limit, and 
the explicit use of CFT constructions rather than direct 
"classification" of spaces of polynomials in the pattern of 
zeroes approach. 



D. Moore-Read construction 

In the MR construction [7|, trial ground state wave- 
functions are constructed as correlation functions in a 
two-dimensional CFT. More precisely, the functions are 
conformal blocks, or chiral parts of correlation functions, 
which depend only on z, not on z, and take the general 
form 

N 

<oio_ w n«(*)i°> ( li4 ) 

i=l 

for ground state, where |0) is the vacuum of the CFT, 
and 

a(z) = e llp{z)/ ^i;(z), (1.15) 

where ip(z) is a free chiral scalar field, and iji(z) 1S a Vira- 
soro primary field with Abelian statistics from some CFT 
( m ore precise requirements will be explained later). 
The coefficient of <p in each exponential is a U(l) charge. 
C_at is an operator of charge —N (so that overall neu- 
trality is satisfied). In MR, On was e" for some 
uniform density p and domain T>, such that it carries to- 
tal charge —N. By including additional fields inside the 
chiral correlator, it is possible to produce wavefunctions 
describing quasiholes as well; we will not require the ex- 
plicit form of these in this paper. 

The operator product expansion (ope), essentially the 
short-distance expansion of correlation functions, plays 
an important role in CFT. The properties of the opera- 
tors a as their coordinates become equal determines the 
algebraic structure of the theory. This information is 
then related to the entanglement subspaces already dis- 
cussed. We will provide very detailed connections be- 
tween the two later in this paper. 

There is a completely general duality, explained in Sec- 
tion |TlT] below, that allows any spaces of symmetric poly- 
nomials that are characterized by "vanishing conditions" 
(or equivalently as the zero-energy spaces for some spe- 
cial Hamiltonian) to be constructed in a similar form 
(again see Ref. [23j], and also the earlier In an 

informal account, we introduce an infinite set of oper- 
ators £_„, n — 0, 1, . . . , all of which commute, and 
let £(z) = Sn>o£-™ z "' they form an algebra A. We 
suppose that there is a vacuum |0). Then the vector 



space spanned by the vectors Yli=i £-n< |0) is dual to the 
space of symmetric polynomials in N variables (much as 
if each £_„ is a creation operator for a boson in the nth 
orbital in the LLL). But now we may also impose local, 
z-independent algebraic relations among the £(z)'s, such 
as t;(z) 2 = 0, or ^(z)d^/dz — 0. In general, each rela- 
tion can be a linear combination of products of £ and its 
derivatives all evaluated at z, and must be homogeneous 
both in £ and in derivatives. Clearly, if such a relation is 
multiplied by any number of £( z i) at arbitrary positions, 
or if linear combinations are formed, other valid relations 
are obtained — thus the combinations that are to be set 
to zero form an ideal I in the algebra A with which we 
started. Then the symmetric polynomials that are dual 
to the quotient A/T acting on |0) must obey vanishing 
conditions as some number of z's come together, corre- 
sponding to these relations. Conversely, for any vanishing 
conditions (of the usual translation and rotation invari- 
ant form) determining spaces of symmetric polynomials, 
such an ideal I can be found. So part at least of the struc- 
ture of the CFT construction above is general [with a(z) 
in place of £(z)]. It is less clear whether the additional 
structure that the polynomials form a module over the 
symmetric polynomials and over the positive part of the 
Virasoro algebra (which is automatic in the MR construc- 
tion), can be shown to be required, though as mentioned 
above, the first at least seems to be necessary in order for 
many polynomials satisfying the vanishing conditions to 
exist. In any case, finding a set of relations that generate 
such an ideal I corresponds to finding special Hamilto- 
nians with a minimal set of terms. Then the search for 
special Hamiltonians takes on a truly algebraic form. 

In recent work, it has been emphasized that the func- 
tions produced by the MR construction for a given CFT 
can be viewed using the operator formalism for CFT. 
Then the operators a(z) are chiral vertex operators, act- 
ing in an auxiliary Hilbert space, which is that of the 
two-dimensional Euclidean (i.e. one space, one time di- 
mension) field theory, not the "physical" Hilbert space of 
particles in two space dimensions. Then such a construc- 
tion corresponds to a continuous version [29] of the pop- 
ular idea of a matrix product state (MPS) — continuous, 
that is in the sense that ZiS are continuous variables, and 
a continuum, and in fact chiral, CFT is used. In partic- 
ular, the (infinite-dimensional) Hilbert space of the CFT 
corresponds to the auxiliary space of the MPS construc- 
tion. Closely related ideas are proving fruitful in on-going 
numerical work on the trial states \4U [421 . 



II. ANALYSIS OF A SPECIAL HAMILTONIAN 
WITH A CONTINUOUS PARAMETER 

In Ref. [lH, it was shown that there is a simple form 
for the (one-parameter family of) ground state wavefunc- 
tions in the MR construction applied to the J\f = 1 SCFT 
[where the field ip(z) — G(z), the superconformal cur- 
rent] , which satisfy a corresponding family of "vanishing 
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conditions" as three coordinates become equal. This re- 
sult can be interpreted as giving the ground state wave- 
function for corresponding special Hamiltonians, whose 
form can be explicitly determined. In this section we 
construct explicitly all the zero-energy wavefunctions for 
these Hamiltonians and enumerate them. We empha- 
size we can view this as a direct analysis of the spaces 
of functions determined by the given allowed behavior of 
the three-particle functions that is one of this family of 
forms (or of zero-energy states for a corresponding spe- 
cial Hamiltonian) , with no use of CFT, other than that 
it motivated the parametrization in terms of the central 
charge c. Though the analysis is somewhat similar to 
other examples (such as those in Ref. [U, S3, the 
present example seems to be the simplest in which there is 
a continuously-varying parameter in the functions. The 
section is divided into two parts: the first gives results 
for the plane geometry, the second for the sphere with 
variable number of flux, which includes the plane results 
as a special limiting case. 

A. Zero-energy states in the plane 

We have seen that for three particles, the symmetric 
polynomials of degree six in the internal or relative co- 
ordinates span a two-dimensional space. Then we may 
consider as a Hamiltonian the sum of projection opera- 
tors onto all states of any group of three particles with 
internal angular momentum (for the three) less than six, 
and onto one of the two linearly-independent states of 
internal angular momentum six. Using Ref. (22[, we will 
parametrize the allowed behavior using first the four par- 
ticle (non-symmetric) function 

x(zi,z 2 ; Z3, Z4) — Azf 3 z 24 zf 4 z 23 + z \ 3 z\ 4 z\ 4 z\ 3 , (2.1) 

where = Z{ — Zj for all i, j, and A = c/3 — 1 where c 
and A can be complex numbers. As zi, z%, z 3 — ¥ Z, this 
becomes 

X~ X 3(zi,z 2 ;z 3 )(Z-z 4 ) 6 , (2.2) 

where 

X3(zi,z 2 ;z 3 ) = lim z 4 & x(zi, z 2 ; z 3 , z 4 ) (2.3) 
24—^00 

= Az 13 z 23 + z 13 z 23 . (2-4) 

When symmetrized over 1, 2, 3, this parametrizes the 
allowed (internal) behavior for three particles in relative 
degree 6 (they must vanish in relative degree less than 
six). Apart from overall normalization of the polynomial 
(which involves integration over all space, and depends 
on the geometry used), this covers most of the space of 
quantum states in the two-dimensional space, which are 
uniquely parametrized by a point on the sphere CP 1 . 
Namely, for finite complex c^oowe cover all except one 
pole. By first dividing by c, we can cover all except the 
opposite pole (c = 0) of CP 1 . 



We will first study the behavior of symmetrized func- 
tions when two coordinates coincide. (The exposition will 
be somewhat terse, as the approach is similar to that in 
the earlier physics literature, in Ref. [HI, S3, S3-) For the 
three particle unsymmetrized function X3 , °ne finds that 
for z\ = z 2 it is given by (A+l)zf 3 . This vanishes if c = 0. 
In this case, the analysis of the space of allowed functions 
is straightforward: allowed functions in N > 2 variables 
must vanish when any two coordinates are equal, so must 
be of the form D^f where Dm = W^jizi ~ Zj) 2 is the 
discriminant (the square of the Laughlin- Jastrow factor) 
in N variables, and / is an arbitrary symmetric poly- 
nomial. In this case the allowed (symmetrized) internal 
three-body behavior in degree six is the same as D 3 . The 
functions of this form are the Laughlin ground and quasi- 
hole states for bosons at v — 1/2, as described in his 
original paper. For N = 2, the allowed functions are all 
of the symmetric polynomials in two variables. This dif- 
ference from the functions for the pseudopotential model 
may seem surprising, but it is allowed by the compatibil- 
ity conditions mentioned in Section II A 31 

From here on, we restrict to c 7^ (with c = 00 handled 
as described above). Then by referring to Ref. [Hj], we 
expect the ground state (the lowest degree polynomial 
of zero energy) to have filling factor v = 1/3. This is 
determined by the degree of the ground state in each 
Zi, and a number called the shift S, which is 6, that is 

A = 37V - 6. (2.5) 

Then lirnjv-»oo N/N^ = v. For c ^ 0, the functions do 
not generally vanish when two particle coincide. We con- 
sider the map G\ which sets two coordinates (without 
loss of generality, the last two) equal to Z\. It maps a 
polynomial from the space of symmetric polynomials 
in N variables into the space Ajv_2 <8> Ai of symmetric 
polynomials in Z\, . . . , Zn—2 times a symmetric polyno- 
mial in Z\. This map has a kernel, that is, the space 
of functions it annihilates. The functions in the kernel 
kerCi are again of the form D^f, for / an arbitrary 
symmetric polynomial (we write -D/vAjv for the space of 
such functions, and similarly in other examples). More 
generally, consider the maps C' m which set the last 2m 
coordinates equal in pairs, to Z\, . . . , Z m . That is, C m : 
f(zi, ...,z N )-> f(zi , z 2l . . . , Z\, Z\ : Z 2l Z 2l . . . , Z m , Zm) 
and maps An —> Ajv_2 TO ®A m . Each of these has a kernel 
in Ajv, and kerC m _i C kerC m . 

We wish to describe the space 7jv (which is a module 
over the algebra A of symmetric polynomials) that is an- 
nihilated by our special Hamiltonian for arbitrary fixed 
c. That is, functions in Jjv vanish in a prescribed way 
in degree six as three particles come together, or faster. 
We define F m = kerC m n I N , for m = 0, 1, . . . , [N/2\ 
(for m — 0, Co — id, with kernel {0}). This immediately 
gives a filtration (sequence of inclusions of subspaces into 
one another) 

Fq = {0} C Fi C F 2 C ■ ■ • -F\Ar/2j QIn = Fin/2\+i- 

(2.6) 
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If we consider c = 0, then F m — In for m > 0. On 
the other hand, for c ^ 0, all the inclusions are strict, as 
can be seen because \3 1S nonvanishing in this limit (or 
maybe better seen below). 

The first kernel is F\ — Jn H DnAn- But as we saw 
above, the polynomials in In have three-particle relative 
behavior in degree six that is not of the form D% unless 
c = 0. For c 7^ then, polynomials in the kernel F\ have 
no term of degree six (or lower) in the differences among 
any three coordinates, and so must be of the form Dn 
times a symmetric polynomial that vanishes when any 
three coordinates are equal. But the symmetric polyno- 
mials that vanish when any three coordinates are equal 
defines a module I N lK , that is the zero-energy states for 
the MR state for bosons at v = 1 [TJ [l7| ■ That means 
we can explicitly write down all functions in F\. 

This generalizes to give a description of F m+ i, modulo 
F m , for all to, including to = [N/2\ if we let F\_n/2}+i = 
In- First, we easily see from the allowed behavior in In 
that, for each to, the image of any function in In (in 
particular, those in F m+ \) under C m is divisible by 

n n^- z *) 6 n (^-^> 12 - ^ 

i<N-2m k<m l<k<l<m 

Hence, now defining the induced map on the quotient 
space C m |F m+1 /kerC m : F m+1 /F m -> A N -2m ® A m , 
which is injective by definition, then we find by a simi- 
lar argument as the one just given (which was the case 
to = 0) that its image is contained in the space of poly- 
nomials of the form 

n ( Zi - Zj ) 2 - n Y[{ Zl -z k f- 

i<j<N-2m i<N-2m k<m 

x J] (Zu-Zt) 12 ■I™* 2m ®A m (2.8) 

l<fe</<m 

for to = 0, 1, . . . , \_N/2\. (Here dots separate products 
of distinct scope, and we reuse indices in distinct prod- 
ucts.) The functions in these spaces correspond to the 
"residues" in earlier papers [THEH. From this, it is easy 
to count the number of such states, as the dimensions for 
the modules /]v-2m are known. This is an upper bound, 
as we have not yet shown that a state exists giving each 
linearly-independent residue (that is, that the maps are 
surjective). It is clear that if they do exist, they are lin- 
early independent. 

To complete the argument first, we will construct one 
function that maps to each of the above residues. First 
(following Ref. 43|), we introduce a labeling of the TV 
particles. The N particles are partitioned into m a > 
(q = 1, 2, 3) clusters of sizes r a , where r a = 2 (a = 1, 
2), r a = 1 (a = 3), so N = 2(toi + m 2 ) + m 3 . One 
particle is assigned to each box in a Ferrers- Young dia- 
gram of at most two columns that corresponds to this 
partition. Then the particle coordinates will be written 

as zjj , where i, j label rows and columns respectively 
(as for a matrix i increases down the rows, while j in- 
creases to the right along the rows) of the rectangular 



block consisting of m a rows of length r a . That is, the 
first mi rows form the first block, the following m 2 the 
second, and the remaining to 3 rows of length 1 form the 
third block. Thus the rows of the diagram are ordered 
with (a,i) above (pt',i'), written (a,i) > (o/,i'), if either 
a < a', or a = a', i < i'. We will deal with the MR 
part at the same time as the other parts. It will turn 
out that toi = to in the above residues, while TO2, TO3 
play a similar role for the MR part. The assignment of 
particle coordinates to will be summed over in the 
end, restoring the symmctrization. 

The wavefunctions for given m a and sets of integers 
rij > (1 = 1, . . . , m a , and a = 1, 2, 3) are 

S z \ II *& ) ,& ) i#U 1 }> (2-9) 

- n n n^M^s )- 

a'=2,3 (l,i)>(a',i') 3=1 
(l,mi)>(a,i)>(a',i') j=l 

x^- 2mi ({ Z r^^^ > 

Here j + 1 is treated cyclically mod r a , so zf^ +1 = 

z^i, DN-2m 1 ({ z ij'~ 2 ' 3 ' > }) i s the discriminant in all the 
variables of types a — 2 and 3, the elementary sym- 
metric polynomials e;({^j° 1 }) are in the variables 

Sj=i one ^ or eacri row °f tyP e a - To see that these 
belong to In, first notice that when two coordinates are 
equal, the polynomial under the symmctrizcr S z vanishes 
unless they both have a = 1 and i is the same, that is 
they are members of the same pair, which we term an 
"unbroken pair" when a = 1. When three coordinates 
approach the same value, the function either tends to 
zero as in these variables, which is the allowed form 
prescribed, or else vanishes faster. If coordinates are set 
equal in pairs until there are to pairs and doing so again 
would make the polynomial vanish, then all these pairs 
are unbroken (a = 1 for each) and to = m±. Further, the 
resulting function (consisting of "half-broken" pairs with 
a = 2, and "unpaired" particles with a = 3) has the form 
of (|2.8[) . in which a polynomial in c an be iden- 

tified using similar results from e.g. Ref. [43|. This com- 
pletes the argument, and as a whole it shows that these 
functions form a complete and linearly- independent set. 
The ground (lowest degree) state, which corresponds to 
mi = N/2, to 2 = TO3 = 0, was found in Ref. [22|], and 
has degree N$ = 37V — 6, as mentioned above. 

The number of states for each degree (angular momen- 
tum) can now be derived either from the explicit basis, 
or from the residues, either for the whole thing or using 
results for J™ m] . The degree of each function, omitting 
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the symmetric polynomial factors, is 
6mi(mi — 1) + 6mi(2m 2 + m 3 ) + 2m 2 (m 2 — 1) 
+ (2m 2 + m 3 )(2m 2 + m 3 — 1) + 2m 2 m 3 + m 3 (m 3 - 1) 



-7V(7V - 2) + 2m 2 + m 3 + 



1 



(2.10) 



2 v / 2 

where we eliminated mi in favor of N = 2mi + 2m 2 + 
?7i 3 . Then the character or generating function for the 
dimensions of polynomials at each degree in In is 



ch„ I N = 



(2.11) 



\N(N-2) 



E 



2 ,rn 3 >0:2 m2 +m 3 < 
(-1)»»3 = (-1)" 



^2m 2 + 5m 3 (iTi3+2) 
--2m 2 -m 3 (g)m 2 (</)r 



where (q) m was defined in Eq. (jl.lOp . This gives the 
complete counting of zero-energy states for our special 
Hamiltonian, for any finite number of particles in the 
plane. 

As a check on ch q In, we can calculate the character of 
I 3 , which should be all symmetric polynomials, except for 
those of degrees 0, 2, 3, 4, 5, and one at 6 in the relative 
variables (times center of mass factors). The formula 
for ch/ 3 does agree with this; it can also be written as 



q 5 )/(q) 3 . Also ch q I 2 = 1/(9)2, 



ch/ 3 = q 6 (l + q + q- 
as it should be. 

We can take the limit N — > 00 fixing attention on the 
zero-energy states of fixed degree relative to the ground 
state simply by taking the limit of this formula with m 2 , 
ro 3 fixed. That is 



N 



lim q-i N ( N -Vch q I N 



E 



(q)oo{q)m 2 (q) m 3 



(2.12) 

u; 1 1 (/ i ,,, !(/),,, 

m 2 ,m 3 >0:(-l) m 3 = (-l)J' 

The sum is still restricted in the parity of to 3 , which 
depends on whether N — > 00 through even or odd val- 
ues. This result has factored into separate factors for 
each m a . The factor l/(g)oo = hm m ^ oc l/(g) m is the 
limit of the character for the symmetric polynomials, any 
one of which may be multiplied into any function in the 
space, and in FQHE often represent edge excitations in 
the charge sector. The remainder of this expression will 
be interpreted in terms of CFT later. 

An extension of what we did here would be to consider 
the states with some quasiholes at the center of the disk. 
For finite size these are merely special cases of the full 
set we obtained above. But when taking N — » 00 we can 
obtain different limits, and then we get the characters for 
other sectors of the theory. In this way we can determine 
all the sectors without considering either the sphere or 
torus. 



B. Location at a phase boundary 

In the early stages of this analysis, we saw for c ^ 
that if we require the functions to vanish when any two 



coordinates coincide, then there is no possible allowed be- 
havior in relative degree six (or less) for three particles. 
It followed that the allowed functions are of the form Dn 
times a function in /^ R . These are the zero-energy func- 
tions associated with the MR state at v = 1/3, and are 
independent of c. Hence, a special Hamiltonian with a 
two-body projection operator (with positive coefficient) 
in addition to the six three-body ones (including the sin- 
gle c-dependent one in degree six) already included pro- 
duces for its zero-energy spaces the functions in DnIn 1R - 
The fact that there is a continuous family of Hamiltoni- 
ans that produce v = 1/3 MR states for bosons does not 
seem to have been emphasized previously to our knowl- 
edge. 

We will presume that the latter Hamiltonians generally 
have a gap in the bulk energy spectrum, and so lie in the 
MR phase (for any c / 0). But when the coefficient of 
the two-body term becomes zero, the zero-energy spaces 
are much larger, and the energy spectrum of bulk exci- 
tations should become gapless in the infinite-size limit, 
because for generic c this system corresponds to ampli- 
tudes from a non-unitary CFT. We conclude that, for 
any value c ^ 0, there is a phase transition when the 
coefficient of the two-body term is reduced to zero. The 
transition appears to have the same behavior for all non- 
zero c, as the zero-energy states have the same multiplic- 
ities, however their wavefunctions do depend on c. These 
states can be interpreted as containing varying numbers 
of the "unbroken pairs" which are more tightly bound 
than the "half-broken pairs" present in the MR ground 
state. When the coefficient of the two-body term be- 
comes negative, the system presumably enters a phase 
different from MR, the nature of which is not obvious. 
However, we may speculate that its ground state con- 
sists of the more tightly bound pairs produced by the 
attractive two-body interaction, something like a strong- 
pairing phase. The presence of the phase transition is 
similar to that in some other examples discovered earlier 

[III EE El. 



C. Zero-energy states on the sphere 

It is simple to repeat the earlier analysis of the residues 
(z\, Z2 — > Z\, etc) for c ^ for the case of the sphere, in 
which the functions of z should be the same degree for 
each particle i = 1, . . . , N. For the residues, the func- 
tions take a Haldane-Laughlin-like form in the positions 
of clusters z[ a \ i = 1, . . . , m a for a = 1, 2, 3 (the zf^ 
are simply the coordinates of the unpaired particles) and 



quasiholes w 

» 



(a) 



I = 1, 



n a for a = 1, 2, 3. Here 



the w)" 1 appear in factors Y\ i i{z[ a ^ — w]" 1 ') for each a. 
When counting the degree of the wavefunction in each 
particle variable z%, one must realize that each z\ a ^ for 
a = 1, 2 can come from either of two particles, so con- 
tributes only a half to the total degree in Zi. Then one 
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finds that the wavefunctions are of the same degree 



range of I depends on a, j: 



provided 



= {6N + n)/2 


-6 (2.13) 


I = 2m 3 + l, 


. , 2m 3 + n 3 - 


Y2 (a 


= i,i 


= 1), 






I = 2m 3 + n 3 + 3, 


m 


(a 


= i,i 


= 2), 


is even) in each Zi 


regardless of type a 


/ = 1, 


n 2 /2 


(a 


= 2,j 


= 1), 


I = na/2 + 1, 


n 2 


(a 


= 2,j 


= 2), 






1 = 1, 


■ , n 3 


(a 


= 3,j 


= 1). 


ri\ = n, 


(2.14) 










(2.19) 



™2 



n 
1 



"3 = ^"-a 



m 3 . 



(2.15) 
(2.16) 



The last line follows from n\ — 2(m 3 + n 3 ) + 4, and 
agrees with a corresponding relation in the MR state in 
the formulation of Ref . [lij ■ 

To write down the complete set of states such that one 
produces each residue, let us first note that the factors \ 
can be used to couple a pair of type 1 to a pair of type 
2, instead of two factors \3 as used above. There is no 
problem viewing these as functions on the sphere. The 
factors X3( z i) 7 - 2 ; z 3 ) are still required to couple type 1 to 
type 3 (unpaired particles). They appear problematic as 
the two terms are not the same degree in z\ and z 2 , and 
are of a different degree in z 3 . The first problem can be 
fixed by incorporating two quasihole factors directly into 

X3- 

X3(zi,z 2 ;z 3 ;w 1 ,w 2 ) = 

A{z x - z 3 f{z 2 - z 3 f{ Zl - Wl )(z 2 - w 2 ) (2.17) 
+ (zi ~ z 3 ) i (z 2 - z 3 ) 2 (z 2 - w 1 )(z 2 - w 2 ). 

As Z\, z 2 — > Z, this tends to the \3 above for the plane, 
times (Z — W\)(Z — w 2 ). The lack of symmetry between 
Wi and w 2 is not important, because of the eventual sym- 
metrization over z^s. The relation n\ = 2(m 3 + n 3 ) + 4 
shows that there will always be more than enough quasi- 
holes of type 1 to allow this to be done for all \3 factors 
that appear. \3 ca n be multiplied by additional factors 
Wl—i 2 {z\ — u>i)(z 2 — Wi) to obtain the same degree in z±, 
z 2 , as in z 3 , fixing the second problem also. 
The wavefunctions on the sphere are 

[a=l,2 (l,t)>(a,i') 

n v / (1) Jl). z (3). (1) . 
X3\ z il ) z i2 ! Z i'l i W 2i>-n W 2 i> J 

(l,t)>(3,t') 

n n^-^si-!-^)- 

(l,mi)>(a,t)>(a',i') j=l 

x^_ 2mi ({,r 2 ' 3) }) nnn (aj) ^ - ^ ah 

a i,j I 

The notation is identical to that for the plane, except for 
the introduction of the X3 above, and the final product 
of quasihole factors, in which Yli appears, where the 



Each of these functions produces just one of the residues 
when particles coincide in pairs. As the residues are sym- 
metric under permutations of w ; (for fixed a), it is per- 
missible to symmetrize over the w\°^ of each a in the full 
wavefunctions . 

Now we calculate the character, which here is the 
group-theoretical character for the representation of the 
rotation group of the sphere acting on the zero-energy 
states, or in other words the generating function that 
counts the number of states of each L z . We use the q- 
binomial coefficient in the form 



m a + n a \ _ (m a + n a ) q \ 
{m a ) q \(n a ) q V' 



(2.20) 



where the q- factorial [n) q \ = [n] q [n — l] q ■ ■ ■ [l] q , and the 
g-deformed integers are 



-(n-l)/2 + -(n-3)/2 _^ h Jn-l)/2 



q -n/2 _ q n/2 
q -l/2 _ q l/2 ■ 



(2.21) 



(This differs from what was used in Ref. [44j, but is 
more convenient for rotationally-invariant situations.) 
For q 1, the g-binomial reduces to the usual binomial 
coefficient. Then by counting the symmetric polynomi- 
als in the that occur in the residues, the generating 
function is 



tr q 



e n 

mi,m2,m,3:2mi+2m2+m3— N a— 1,2,3 

(2.22) 

The results for counting states in the infinite plane 
can be recovered by letting n — > oo, as follows. In the 
plane, by stereographic projection from the sphere, — L z 
corresponds to angular momentum in the plane (total 
degree), but we must shift — L z by adding — 
N(6N+n) /4—3N. Combined with the highest L z (which 
occurs for N even, mi = N/2, m 2 = m 3 = 0), which is 
L z = Nn/4, we recover the degree of the ground state in 
the plane, ^N(N — 2) (independent of n). After extract- 
ing this factor, the g-binomials take on the definition used 
in Ref. Q, in which [m] q -> (m) q = (1 - q m )/(l - q), 
and this is convenient for taking limits such as n — > oo at 
fixed N . In this way the generating function (character) 
in the plane, Eq. (|2.12[) . can be recovered. 

Returning to the sphere, we can view the g-binomial 
in mi, ni as describing the "positional" degeneracy of 
what would be the quasiholes if the bulk energy spectrum 
were gapped, as in a topological phase: there are n = n\ 
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quasiholes of charge 1/6. The other factors describe the 
degeneracy for fixed positions. This depends on 777,2, 7773, 
while mi = (2m2 + m + 3)/2 is the number of unbroken 
pairs. The total degeneracy for fixed positions is 



fm 2 + n 2 \ fm 3 + n 3 

m2,m3 ^ 2 / \ 



n 3 



(2.23) 



where the sum is over 777,2 and 7713 such that N—2rri2 — 7773 
is even and > 0, and 77 2 , 773 are determined by n and 7713 
as above. In a sensible phase, this would approach a finite 
limit as N (or AT^) — » 00 with 77 fixed. In particular, 



7773 + rt 3 
™3 



n 2 /2 
7713 



(2.24) 



corresponds to 777,3 fermions that can occupy 772/2 states. 
This factor is never large, and 777,3 must be less than or 
equal to n 2 /2 = n/2 — 2. But the binomial coefficient 
in 7772, n 2 corresponds to 7712 bosons in 772 + 1 orbitals, 
and as N — > 00 with 77 2 > there is no upper limit 
on m 2 . That is, the number of half-broken pairs can be 
arbitrarily large as N — > 00. Hence the total degeneracy 
for fixed positions of the quasiholes becomes arbitrarily 
large, similar to the permanent and Haffnian cases 
[45I ] . This behavior cannot occur in a gapped (topological) 
phase of matter. 



III. AMPLITUDES AND MODULES 
A. CFT — operator methods 

In this section, we first set up definitions and basic 
general arguments that will be used in the following. The 
starting point is the general form of the MR construction 
of trial wavefunctions, eq. (|1.14j) for the ground state. 
Then we describe our first main result. 

In this paper, the main examples are paired states for 
which the operator product expansion (ope) of ip reads 



1 



+ 



(3.1) 



as z — > 0, where is the conformal weight of ip, and . . . 
denotes lower-order terms (not only non-singular terms) . 
In our main examples, either (i) tp is the field at the 
corner of the minimal block in a BPZ minimal model 
M(p,p'), that is ip = </>(i iP ) or 4>( p >.i) m the p' mini- 
mal model (we use the conventions of Ref. Q), or (ii) ip 
is the superconformal current G(z) in an J\f = 1 super- 
conformal field theory. The first of these includes, but 
is more general than, the </>(2,i) field used by FJM. (For 
the Laughlin states, the analysis is similar but easier.) 
We will, however, extend the general results to the case 
of generalized parafermionic CFTs, or clustered FQHE 
states, of order k = 2, 3, 4 . . . ; in these CFTs, when 
the ope of ip's is iterated, k ip's are needed in order to 
obtain the identity operator on the right hand side, as 



for the k = 2 case above (the case k = 1, where ip = I, 
the identity operator, gives the Laughlin states only). In 
all cases, the right hand side of the ope of ip with itself 
is assumed to contain only operators of the same scaling 
dimension as the operator in leading term, plus integers, 
and similarly for the iterated opes; ip should be a "simple 
current" (these conditions ensure that, for suitable values 
of h and h^,, the wavefunctions we obtain are indeed just 
symmetric polynomials) . The assumed behavior of the it- 
erated ope's of ip leads to restrictions on the value of h$ 
(we will not need the details). For paired (k = 2) cases, 
single- valuedness and symmetry implies that v — 2h^ 
must be a non-negative even integer. In particular, one 
can consider v = l/(2h^,) so the wavefunctions are non- 
vanishing when two coordinates coincide. For general fc, 
we require that the conformal weight h = v~ l /2 + of 
a, be half a positive integer so that a and its adjoint gen- 
erate a chiral algebra, or vertex operator algebra (VOA) 
!). 50j, as we will discuss below. In some cases 2h is 
odd, even though the wavefunctions that are produced 
describe bosons — the chiral algebra generated by a and 
a* violates the spin-statistics theorem (this occurs only 
when ip is part of a non-unitary theory, such as the case 
of the Gaffnian for bosons) . We can also obtain antisym- 
metric wavefunctions for fermions in a similar way, either 
by relaxing these conditions (leading to chiral superalge- 
bras, but still obeying or not obeying the spin-statistics 
relation), or simply by multiplying the wavefunctions for 
bosons by the Vandermonde determinant Ili^^i — z j)- 
We have the charge current operator j(z) — idip, with 
operator product 



mip(O)/ \/v 



intfi(0)/y/u 



(3.2) 



Thus j measures the charge, up to the factor ^/v. Being 
a CFT, the theory has a stress tensor 



T(z) = T v (z)+T i ,(z), 



(3.3) 



which is a sum of the stress tensors T v and of the 
theories to which ip and ip belong. We have T v (z) — 
— i(c/V) 2 . Finally, we have a vacuum |0) of the combined 
CFT, which is the tensor product of the vacuua of the 
ip and ip theories, |0) = 10^,) ® |0^,). We emphasize that 
in this section, all states and operators are those of the 
two- (1 + 1-) dimensional CFT, not the 2 + 1 dimensional 
many-particle system. 

The mode expansion of a can be defined by 



a(z) = ^ a n z 

neZ 



(3.4) 



as z — > (acting to the right). The modes can therefore 
be extracted as 



a„ = — !— d> dz z n 1+h a(z). 
2-ni r 



(3.5) 



The integral is over a circle of radius r. As h is an integer, 
and there is no singularity at the origin if no non-trivial 
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operator is inserted there, we must have 

a„|0) = (3.6) 

for n > —h, and we can identify the state a(0)|0) = 
| a) = a~h\0). All products of z-dependent operators are 
radially-ordered, unless otherwise stated. To obtain the 
product of modes in any given order, one uses the corre- 
sponding ordering of the radii of the circles in the integral 
defining each mode. 

The mode expansion takes the same form for any con- 
formal (Virasoro primary) field, with h replaced by the 
conformal weight, and the mode indices n must obey 
n + h £ Z. In particular, the modes j n of the weight- 
1 U(l) current j(z), and L n of the weight-2 stress tensor 
T(z), obey the usual relations, in particular the Virasoro 
algebra relations for the modes L n — L nip + L n ^ of T in 
radial quantization Also, the modes obey relations 



with a(z), 

\j n ,a(z)] = ^ r z n a(z), (3.7) 

[L n ,a(z)} = h(n+l)z n a(z) + z n+1 ^^. (3.8) 

oz 

Finally, we also have 

j„|0) = (n > 0), (3.9) 

L n |0) = (n > -1), (3.10) 



and we refer to the sets of modes in the last relations 
as the positive half of the current and Virasoro algebras, 
respectively. 

For our main examples, in which k = 2, it is clear that 
the conjugate field of a (which creates an "antiparticle" 
from the vacuum, if a creates a "particle") is 

a*{z) = e- l{p{z)/ ^i>(z). (3.11) 

This construction works because tp is its own conjugate, 
a* has the same conformal weight h as a. Clearly a* can 
be obtained by the local operation of ordinary complex 
conjugation, with f(z), ip(z) treated as invariant (self- 
conjugate). 

We now assert that this operation can be viewed as the 
radial adjoint (here we differ slightly from some standard 
treatments, such as Ref. @): 

a*(z) = a(z) t . (3.12) 

For any conformal field, the radial adjoint is defined using 
radial quantization, and the radial coordinate plays the 
role of imaginary time (because of the Euclidean signa- 
ture on the two-dimensional spacetime; note also that the 
definition of the adjoint operation depends on a choice of 
origin). That is, we first map conformally to the cylin- 
der, with the axis of the cylinder as imaginary time; the 
map is z i — y — ilnz. We then apply the adjoint opera- 
tion, which strictly speaking is only the usual operator 
adjoint at \z\ = 1 (imaginary time equal to zero); for 



non-zero time, it involves reversing the sign of imagi- 
nary time (in addition to complex conjugation), as in 
the imaginary-time formalism in many-body theory [48[ . 
That is, — zlnz is invariant in this step. Finally, we ap- 
ply the inverse conformal map back to the plane. The 
conformal map to the cylinder means multiplying by a 
factor z h , and by the inverse for the inverse map. Thus 
in this point of view, = z -1 , except in the z~ h factors 
which are invariant. Putting all these steps together, we 
have 

a(z)t =£ a t*»-\ (3.13) 

n 

In modes, using the standard form of mode expansion for 
a* or a\ this then implies 

( a t)„ = ( _ n )t. (3.14) 

(This last formula is standard @.) The latter formulas 
apply to all conformal fields, for example ip(z), and also 
to j(z), T(z) which are both self-conjugate, so that j- n = 
jj, and so on. 

Usually adjoints are defined relative to an inner prod- 
uct, but in CFT we frequently do the reverse. We take 
the preceding formulas for the adjoint operators, together 
with the inner product of the vacuum (0|0) = 1. Inner 
products of states defined as any finite collection of modes 
of the fields acting on |0) can now be defined using the 
preceding formulas, and using the commutation relations 
to simplify expectation values of the modes. This deter- 
mines the inner product on the space of states, which is 
Hermitian but not necessarily positive definite or non- 
degenerate. 

We define 

\0 N ) =e ^(o)/^| ), (3.15) 

or in an obvious notation |0jv) = \0 v n) ® |0^). Also we 
can define an out vacuum of charge —N by 

(0 N \ = lim ^v>( | e -^(z)/v^ ( 3 . 16 ) 

Z— >oo 

[We can do similarly for out states of other conformal 
fields, using twice their conformal weight in place of 
N 2 /v, thus (a\ = lini2_ i . 00 z 2h (0\a* (z) .] Then we have 

(0 N \0 N ,) = 6 NjN ,. (3.17) 

We can now define the chiral algebra, or vertex op- 
erator algebra (VOA) [47| . of our theory as the algebra 
generated by the modes of fields generated by repeated 
ope's of a(z) and a*(z). (We do not mean to imply that 
a chiral algebra and a VOA are the same concepts when 
defined rigorously, however for our purposes greater pre- 
cision will not be required.) We will argue that it includes 
the current operator j(z) and both stress tensor opera- 
tors, T v {a) and T^{z). 

For the ip theory, we define the space of states Tn 
(N € Z), to be the Fock spaces generated from |0 v at) by 
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the action of the modes j n (n < 0). For the ip theory, 
we define the space M (Mi) to be the space of states 
generated by the action of an even (odd) number of ^;'s 
on the vacuum |0^). In the Virasoro minimal model case, 
Mo is the vacuum (1, 1) irreducible Virasoro module, and 
Mi is the irreducible Virasoro module generated by 
In the SCFT case, we will discuss the structure of the 
spaces later. Then, as in Ref. |23[, the space V of all 
states spanned by a and a* acting on |0) is V = (§) N Vn, 
where 



V N = 



T N ® Mo (N even) 
F N <E> Mi (N odd). 



(3.18) 



The inclusion V C W Vn is clear, but the equality re- 
quires proof. First, we notice that V is a module over (i.e. 
a representation of) the VOA. By the state-operator cor- 
respondence (or the theory of VOAs) d H, V "has 
the structure of" the VOA, in particular there is one 
(linearly-independent) state for each non-zero, linearly- 
independent operator in the VOA, and the Virasoro level 
and the charge of the states and operators agree. 

It is not too difficult to see that V as defined here is 
indeed generated by a(z) and a(z)* . Essentially, from 
the ope's one can extract j(z), and hence separate out 
the spaces J-n- The generalization from order k — 2 to 
k > 2 is straightforward, involving k modules Mo, . . . , 
Mfe_i in place of the two above. A more detailed proof 
for all cases is given in Appendix |A"1 

The mode operators a n commute, because the oper- 
ator product a(zi)a(z2) is non-singular. The operators 
a generate a commutative subalgebra U (the principal 
subalgebra) of V, and the subspace spanned by these op- 
erators acting on |0jv) will be denoted V N : 



V N = U\0 N ) 



(3.19) 



Then we have a chain of embcddings (injective maps) 



V 4 



V 1 



v° 



v 



(3.20) 



all of which are subspaces of V. The embeddings arise 
because V N ~ 2 includes a(zx)a(zi)\QN-2), and the leading 
part as Z\, z% — > gives |0jv), and hence all of V N appears 
inside V N ~ 2 . We call W — V° the principal subspace (it 
corresponds to the principal subalgebra). There is also 
[23j an operator that we will call Q that shifts the charge 
by 2 (it is part of the operator construction of e 2%v ^ 
and acts as 



Qa n Q~ 



(3.21) 



on the modes. The map Q : V N —> V N+2 is an isomor- 
phism. Q will be used later. The appearances of the 
number 2 in these statements should be replaced by k in 
general. 

It was shown in Ref. [40( and FJM, for the particular 
set of models they consider, that V is isomorphic to the 
direct (or inductive) limit lim > V~ N of the above chain 



of embeddings, defined by 



limV 



-JY 



lim V 



-N 



N 



(3.22) 



JV 



Their proof uses the irreducibility of Vn (Vn is the sub- 
space of V of charge V) under the action of the algebra 
generated by the current j(z) and the stress tensor T(z). 
This holds in our minimal model constructions, of which 
the models in FJM are special cases. A more general ver- 
sion, for all k — 2 constructions, can be proved using the 
algebra generated by ip itself, as can the generalization 
to all k > 0. That is, we have 

Theorem 1: In notation as above, for the order k case, 



V = lim V 



— JV 



(3.23) 



Equivalent ly, V is generated as a VOA by a(z) and 



(3.24) 



(The latter is needed to reach the negatively charged vac- 
uua |O_jv)0 We again defer the detailed proof, which is 
not difficult, to Appendix [A] This is the first main re- 
sult of this Section. (For the Laughlin or k = 1 case, 
the Theorem is trivial, as V is generated by a(z) and 
a(z)* = e ~^(2)/V^ by definition.) 

The principal subspace W contains vectors constructed 
from the commuting operators a(z). W is a module (rep- 
resentation) of the principal subalgebra (i.e. itself, under 
the state-operator correspondence). From the relations 
(|3.7p , (|3.8[) , we have commutation relations for the modes 
of a(z), 



[L ni a m ] = [(h - l)n - m]a m+n , 



(3.25) 
(3.26) 



and see also eqns. (|3.9[) . (|3.10p . W is thus also a mod- 
ule over these positive halves of the current and Virasoro 
algebras. This positive half of the Virasoro algebra in- 
cludes the generators Lq, L±i of the Mobius group, but 
the central charge never shows up. Note that L_i gener- 
ates translations. 

We can describe the structure of W by considering an 
algebra A spanned by products of the commuting inde- 
terminates £_ n , n = 0, 1, . . . (which correspond to a_ n _^ 
by dropping the shift by h). The operators a in general 
obey further relations, and these relations generate an 
ideal I c A. Then we have the isomorphism 



W A/1. 



(3.27) 



We emphasize that the ideal X of relations in A can be 
thought of as generated by local relations among the a's, 
because of locality in the CFT. Examples of such rela- 
tions are a(z) 2 = in the Laughlin v = 1/2 example, and 
a(z) 3 — in the MR v = 1 example, in which ip is a Ma- 
jorana fermion field. The relations must be independent 
of position z (translation invariance), respect the degree 
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d = n i °f the products of £_nS (rotation invariance), 
and respect the charge N (charge conservation). Each 
such subspace Ad,N of degree d and charge N is finite di- 
mensional. Any such expression written as modes of a(z) 
acting on |0) can be viewed as built from a's and deriva- 
tives of a's acting at z = 0. The translation invariance 
of the VOA allows such a relation to be moved to any 
point z. Moreover, there is a natural action of the (pos- 
itive) current and Virasoro modes on A, determined by 
the same relations (I3.25[) , (|3.26p . The subspaces X N (of 
vectors with charge N) must be mapped into themselves 
by these modes. 



B. Dual spaces of polynomials 

We now turn to the "functional model" of W [H, lipj 
which corresponds to wavefunctions for bosons in the QH 
situation. We recall that the dual B* of a vector space 
B is the space of linear maps from B to the complex 
numbers, and the action of a map d € B* on a vector 
b € B can be written as the "dual pairing" (d, b) (no 
complex conjugation is involved). The dual vector space 
of A can be represented by symmetric polynomials; we 
write Ajv for the symmetric polynomials in N variables 
(which are here usually written z\, . . . , zn), and Ad,N 
for the subspace of polynomials of degree d. We have the 
natural dual pairing [23[ 



(f(zi,...,z N ),£- ni •••£-„„) = 
Res Zl= ... zw . =0 /(zi, ■ • ■ , z N )z^ ni ~ 1 ■ 



(3.28) 



-njv-l 



where / is a symmetric polynomial in Ajy , and the residue 
is the usual one of complex analysis, applied to all Zj. 
That is, we simply pick out the part of / (expanded in the 
"monomial basis" of symmetrized products of powers of 
zi) in which the degrees in the N distinct variables are (in 
any order) m, . . . , tin- The dual space W N of Wn (the 
subspace of W spanned by N a's acting on |0)) is then 
the subspace In of A at of polynomials that annihilate the 
module In with respect to the above pairing, 



I N = {f€A N \(f,0=0^€TN}. 



(3.29) 



This means that In can be viewed as the space of sym- 
metric polynomials that satisfy certain "vanishing condi- 
tions" on their behavior as several z^s come to the came 
value; these conditions are derived from the local rela- 
tions in I. 

Equivalently, we may think of In as the space of all 
symmetric polynomials constructed as 



f v (zi, . . . ,z N ) = (v,a(zi) ■ ■ ■ a(z N )\0)) 



(3.30) 



where v is an element of the dual space V* of V, and 
( , ) is the dual pairing with V. Clearly such functions 
inherit vanishing properties from the relations satisfied 
by the a(z)s. This form is very close to what we wish 
to use for wavefunctions, but it will be important later 



to realize that the use of the dual space V* here is not 
always the same as using the space of out states. The 
out states can be viewed as elements of the dual space, 
which act on the space V by the natural (Dirac) pairing. 
Thus the out space can be mapped into the dual space; 
however this map may not be injective. We will now 
discuss this. 

If the inner product (Hermitian sesquilinear form) on 
V is non-degenerate, then we can identify the dual space 
V* with V using the inner product (more precisely, we 
have a conjugate-linear map from one to the other). We 
recall that an inner product is degenerate (or singular) 
if there are vectors that have vanishing inner product 
with all other vectors; such vectors are called singular (or 
sometimes, miscalled null) vectors. The singular vectors 
form a subspace, and the space modulo the subspace of 
singular vectors (formed by "setting the singular vectors 
to zero") possesses a non-degenerate inner product. 

Instead of using the dual space V* to obtain the poly- 
nomials, we can use the inner product and obtain poly- 
nomials 



f w (zi, . . . , z N ) = (w\a(zi) ■ ■ ■ a(z N )\0), 



(3.31) 



where w 6 V. We denote by Jn the space of such poly- 
nomials for N a's. Clearly, Jn Q In- If the inner prod- 
uct is non-degenerate, then Jn = In- But if the inner 
product is degenerate, then Jn may not equal In, as we 
have already noticed above. In finite dimensional vector 
spaces, the space and its dual have the same dimension. 
The space In in each total degree d (we write these as 
Id,N, and similarly for JcI,n) has the same dimension as 
Wn at L eigenvalue d + Nh (which we write as Wcl^n)- 
But Jd,N has the same dimension as the quotient space, 
Wd,N modulo singular vectors lying inside it. Jn can be 
identified as the dual of that quotient space. 

The spaces of polynomials In and Jn at which we 
have arrived have further general properties that reflect 
general aspects of our construction. The right action of 
the current modes j n (n > 0) and the Virasoro modes 
L n (n > — 1) on v G V* give [23[ 

fvj„(zi, ...,z N )= —=s n f v (z x , ...,z N ) (3.32) 

fv-L n (zi,---,z N ) = h(n+ l)s n f v (zi,...,z N ) 

+ £ n {h)(z 1 ,...,z N ) (3.33) 



where for n > 



A' 



i=l 

are the sums of powers (so = A), and for n > — 1 



(3.34) 



i=l 



- 1 — . 

dz,' 



(3.35) 



These follow from eqns. (GO]), ([3T5]). and (|3~TU1) . 

The sums of powers s n (with exponents < n < N) are 
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another set of generators for the symmetric polynomials 
An, and so the space In is mapped into itself by mul- 
tiplication by symmetric polynomials — it is a module 
over the algebra of symmetric polynomials. Similarly, Tjv 
is a module for the positive half (actually n > —1) of the 
Virasoro algebra, acting by the £ n 's. That is, s n lN Q In 
and £ n lN Q In- Both sets of operators raise the degree 
of the polynomial (for n > 0), and usually one might 
expect them to be associated with the raising (negative) 
parts of the algebras; the apparent reversal of sign is due 
to the dual relationships in use here. Identical arguments 
apply if the operators are applied (from the right) to the 
out-state vector (w\ to obtain an action on the functions 
in Jpf. If we think of the adjoints j'_„ and L_„ acting on 
the vector \w), then these are indeed raising parts. 

It is useful to realize that, in terms of the interpretation 
of symmetric polynomials as wavef unctions in the LLL, 
£_X is a derivative with respect to the center of mass 
variable, while j\ acts as multiplication by (N times) the 
center of mass. These obey the Heisenberg relation 

[£- U h]=-N, (3.36) 

corresponding to the projective representation of the 
translation group of the plane. 

Turning again to the question of singular vectors, for 
a degenerate inner product in the context of a CFT or 
VOA, the subspace of V (or of any other module obeying 
the adjointness conditions) consisting of all the singular 
vectors must form a submodule of the VOA. Conversely, 
in an indecomposable module of the VOA (i.e. one that is 
not a direct sum of modules), a submodule must consist 
entirely of singular vectors. Finally, the quotient of V by 
the subspace of singular vectors yields a non-degenerate 
VOA on which the inner product is non-degenerate. We 
note that a non-degenerate inner product is a necessary 
condition both for a rational and for a unitary CFT (but 
not sufficient for either). 

In some situations, examples of which we will study in 
this paper, a VOA is first constructed that depends on 
a parameter, such as the Virasoro central charge c. For 
generic values of c, the space V may be non-degenerate. 
But for special "rational" values of the parameter, sin- 
gular vectors may appear @, Q , and then they have to 
form a submodule of the VOA inside V. The central 
point is that this subspace may intersect the principal 
subspace W. In terms of the algebra A, there is then a 
larger ideal X' (over A) that corresponds to the singular 
vectors. That isld'Ci. Then X'/X C A/I = W 
is the singular subspace of W. Then, as we know, Jn is 
the dual of the quotient of these spaces, 

Jn = (W N /(X' N /X N ))* , (3.37) 

or Jn = (An /X' N )* . In other words, the smaller space 
Jn is dual to a space that is modded out by additional 
(independent) relations among the a's, described by I', 
rather than X. That is, the subspace Jn C In is de- 
fined as polynomials that obey additional (independent) 



"vanishing conditions" as several ZiS approach the same 
point, compared with those defining In if Jn In- In 
particular examples, as a continuous parameter such as 
the central charge c is varied, this phenomenon may oc- 
cur at particular "rational" values of the parameter. 

C. Application to entanglement subspaces of 
ground state 

Now we turn to the applications of the preceding ideas 
to QH trial wavefunctions. In this section we translate 
the preceding results, and use some additional ideas to 
derive statements about the entanglement subspaces of 
ground state (or other) wavefunctions obtained from a 
MR construction. 

The starting point is to write the trial ground state 
wavefunction as in eq. (|1.14j) . following MR. In that work 
the operator O-n was a uniform charge distribution, and 
it produces the Gaussian factors e~' z ' / 4 in the wave- 
functions in the plane (in the limit that the background 
charge distribution extends to infinity — note that this vi- 
olates neutrality). These factors are part of the lowest 
Landau level wavefunctions (in the symmetric gauge), 
and are important when calculating quantum-mechanical 
inner products. In this paper we usually do not con- 
sider these inner products, and so the Gaussian may be 
dropped. It is then convenient to use instead the point 
charge operator 

O-n = e - lNlp{z)/ ^ (3.38) 

and place it at Z — > oo. This means that (the polynomial 
part of) the trial ground state wavefunction is given by 
the "conformal block" 

(*!, ...,z N ) = (0;vK*i) • • • a(z N )\0). (3.39) 

More generally, we may then consider trial functions 
that are obtained as the "amplitudes" we have been dis- 
cussing, that is functions in In (or Jn) for the CFT under 
consideration: 

V v (zi,...,z N ) = (v,a(z!- ~a(z N ) |0)). (3.40) 

The spaces of amplitudes In, obtained as v ranges over 
V* (or Jn similarly, using out states) form spaces of nice 
trial wavefunctions. 

The spaces In are dual to Wn- Using the map Q, 
we can map W — V° to V~ N , and similarly for Wn- 
Then Theorem 1 above implies that as N' —> oo, the 
sum of spaces © w In can be identified as dual to V 
itself. That is, the ground state for some N = N' is 
viewed as the vacuum (approximate vacuum, for finite 
N), and other states, that differ from it in charge and de- 
gree by amounts that stay constant as the limit is taken, 
are viewed as excited; then finally take N' —> oo. In 
other words, if we suppose that the space of amplitudes 
can be identified as a space of zero-energy wavefunctions, 
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and the latter can in turn be viewed as a space of edge 
states, we obtain the fact familiar in examples that the 
edge excitations (including those that change the particle 
number) of the ground state of a disk span the (dual of 
the) vacuum representation of the chiral algebra (VOA) 
used in the construction, as N — > oo. (Such a relation 
of edge CFT with that used in the MR construction was 
proposed in MR [7|.) In particular, the characters of the 
spaces agree, provided we subtract the charge (degree) 
of the ground state wavefunction from the total charge 
(degree) before taking the limit — this corresponds to the 
use of the operator Q. This relation is one of those some- 
times referred to as a "bulk-edge correspondence" . 

Now let us consider the spaces Jn of amplitudes de- 
fined using the inner product on the CFT (for the mo- 
ment, it may be the same space as In), and in particular 
view the ground state in this way, as in eq. f|3 .39|) . We 
let N — Na + Nb, and consider Nb of the a(z)s (say, 
those with i = Na + 1, • • • , Na + Nb) as acting to the 
left. This can be done by taking the absolute values of 
their coordinates Zi larger than all those of the ZiS with 
i = 1, .. . , Na- Also, we may expand the former a(z) in 
modes. Thus we arrive at out states of the form 

Nb 

(0N\Ha nk . (3.41) 

fc=i 

Under the map ("adjoint") from V to out states, these 
correspond to 

N B 

]Ja*_ nk \0 N ). (3.42) 

fc=i 

Comparing again with the first main result above, and 
replacing a in the argument by a* (and reversing the 
sign of the charges), we learn that as Nb — > oo, these 
states eventually span V. This means that applying this 
expansion to the ground state, the ground state decom- 
poses into combinations of amplitudes in Jn a , and as 
Nb — > oo, all of Jn a is obtained. 

This decomposition in fact involves the entanglement 
subspaces. Dividing the particles into those with coordi- 
nates inside some disk of radius equal to one, centered at 
the origin, and those outside, is the definition of a real- 
space cut. We may sharpen the preceding argument by 
imagining that we insert a complete set of states in V, 

J=X>X«I (3.43) 

vev 

(for an orthonormal basis set of v) at radius one. Then 
we arrive by the preceding argument at the same con- 
clusion, now for the entanglement subspaces in the limit 
as the size of part B tends to infinity. The dimensions 
of the subspaces in real-space entanglement are the same 
as those in particle partition. This was proved in Refs. 
[2SH3 by a linear-algebra argument. In the present con- 
text, it is nicer to explain it in terms of functions [28j : 



the real-space cut leads to functions in Na (Nb) vari- 
ables restricted to lie in the disk, part A (or its comple- 
ment for B). But these analytic functions can be analyt- 
ically continued to the whole plane, giving the functions 
in the particle partition. Hence we have shown that as 
Nb — > 00, the particle (real-space) entanglement sub- 
spaces of the ground state in a MR construction coincide 
with (are isomorphic to) the spaces Jn a . This is the 
promised result concerning version (i) in the Introduc- 
tion. Though for Na — > 00 this result can be extracted 
from the recent work [2!|, the proof here makes no use 
of the quantum-mechanical inner products on the poly- 
nomials (which required use of a generalized screening 
hypothesis), and consequently applies to a larger class of 
constructions; moreover, the result here holds in finite 
sizes Na- In the orbital partition at fixed Na, the entan- 
glement subspaces are generally smaller than those in the 
other two forms, though it is possible that the dimensions 
become the same as the limit Na — > 00 is taken. 

In the preceding two results, the use of distinct spaces 
Jjv and Jn in each one was necessary in general, as we 
will see momentarily. But when the inner product on V 
is non-degenerate, the two spaces are identical. We then 
obtain a correspondence of the entanglement subspaces 
(as Nb — > 00) with the spaces of amplitudes (or zero- 
energy states, or edge states) for N = Na particles. This 
is a finite-size form of bulk-edge correspondence. 

Now we consider what happens if the inner product 
on V is degenerate, and there are singular vectors in W. 
The space of amplitudes Iff, to be viewed as spaces of 
"nice" functions, for the given CFT in the MR construc- 
tion is a definition. For the entanglement subspaces of 
the ground state, it is the inner product in the CFT that 
is important, and hence the spaces Jn- This can be seen 
in several ways. For example, the use of a complete set 
of states as above presupposes a non-degenerate inner 
product. Or we may note that the same ground state 
function is obtained if in the MR construction one takes 
the quotient of the VOA by all singular vectors; the sin- 
gular vectors do not contribute to correlation functions 
in the vacuum [H, Q . A particular example is the Af = 1 
superconformal model whose wavefunctions were stud- 
ied in Section [TTJ The ground state wavefunction can be 
determined by the superconformal algebra alone [22| , re- 
gardless of whether the singular vectors are modded out. 

This leads to an important consequence for such cases: 
if there are singular vectors in W, the entanglement sub- 
spaces Jn will generally be smaller than the correspond- 
ing space of amplitudes In (or zero energy states, when a 
special Hamiltonian exists), even as N — > 00. This is true 
in particular for the superconformal examples at rational 
central charge. Some of the pseudoenergies (defined as 
minus the logarithm of the Schmidt singular values) in 
the Schmidt decomposition of the ground state must go 
to infinity as c approaches such a value. In this case the 
correspondence with the spaces of amplitudes (or "edge 
states") is lost. However, the correspondence can be re- 
covered by using instead the quotient CFT when finding 
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the amplitudes or zero-energy states. This will however 
involve a change in the special Hamiltonian, as we will 
now explain. 



D. Zero-energy states and special Hamiltonians 

In this Section, we begin to merge the streams of 
thought in this paper, applying the preceding analysis 
of this section to the problem of solving or construct- 
ing special Hamiltonians, using the AT = 1 SCFT model 
of Section |n] as a specific example. In that section, we 
saw that the short-distance ope behavior of the SCFT 
implies that the entanglement subspace for three par- 
ticles in the ground state begins with one polynomial 
in degree six, followed by higher degrees, in the relative 
or internal variables; thus the other degree six, and the 
lower ones are missing. This means that imposing such 
vanishing conditions to be obeyed by the wavefunctions, 
or equivalently using a three-body special Hamiltonian 
that enforces such conditions on zero-energy states of 
the TV-particle system, produces spaces of polynomials 
for each N that must contain the spaces of amplitudes of 
the SCFT, as obtained by the MR construction that we 
have discussed in this section. 

It is not immediately clear whether the spaces of zero- 
energy states are exhausted by the spaces of amplitudes, 
or whether they also contain other linearly-independent 
states. This issue arises in most examples of special 
Hamiltonians, and we cannot always answer it. But in ex- 
amples in which we can obtain all the zero-energy states, 
we can compare the counting of the states with that of 
the spaces of amplitudes. For the latter, we lack detailed 
results for finite N in general. However, it is a standard 
problem in CFT to compute the character of the vac- 
uum module, which as we have seen corresponds to the 
N — > oo limit of the character of the spaces of ampli- 
tudes, and we can then compare these limits. With that 
information at hand, we will be able to show the equality 
of the spaces even for finite N. 

Returning to the system of Section UH we have em- 
phasized that the full characters include the counting of 
states at different charges from the vacuum. But it will 
be sufficient here to examine the limits for a sector of 
fixed charge N' — N, as N' — > oo, because the interest- 
ing part is the ip CFT. The result in eq. (|2.12l) . with 7713 
even, can be viewed as the sector of zero charge excited, 
relative to the vacuum (the other charge-even sectors dif- 
fer only by a power of q, left over when removing the 
power qi N ( N ~ 2 ) for the ground state with N' particles, 
and the charge-odd sectors also involve summing over 7713 
odd instead of even). The expression has factorized into 
1/(9)00, which as we mentioned corresponds to excita- 
tions of the charge sector, and two other factors, each 
of which is a sum. They can be expressed in a different 
form as follows. We use combinatorial identities of Euler 
[4^ . By specializing these general results, we can obtain 



for the TO2 sum 



E 

fc=0 



21, 



1 



(<7)oo 



(3.44) 



For the sum over 7713, we first leave it unrestricted, and 
obtain two identities, 



00 

E 



(±i)V 



kf 2 +f 



ir=i(i±g r - 

l±q 1 / 2 



(3.45) 



We may produce identities for the sums with m 3 even 
(odd) by taking half of the sum (difference) of the ex- 
pressions with the plus and the minus. 

If the simple factors on the right hand side are omit- 
ted (that is, 1 — q in the first, and 1 ± q 1 / 2 in the oth- 
ers), these products have familiar interpretations as the 
partition functions for chiral bosons and Majorana-Weyl 
fermions, respectively. But with those factors, a different 
interpretation is more useful. 

Combining the three factors in the full charge-zero 
character, it will be convenient again to consider the un- 
restricted 777,3 sum, weighted with (ztl)™ 13 (a "Z2 graded" 
character). If we write V for the space of states in the 
limit with the charge sector removed completely (the odd 
and even parts have to be combined carefully with the 
charge sectors, as we have seen), then we have 



ch± V = ch 9i+ i V ± ch g ,_i V 



(i-?)lELi(i±«r-*) 



(«)oo(l±gVa) ■ 
(3.46) 

This character can be identified as that of the Kac mod- 
ule for the vacuum module of the M = 1 superconformal 
algebra (in the Neveu-Schwartz sector) . This is obtained 
as follows. The Verma module [8, 9] consists of all for- 
mal descendants of the identity under the positive modes 
L- n , G_ n+ x/2 of the stress tensor and superconformal 
current, respectively, where 71 > is an integer. Using 
standard arguments (the Poincare-Birkhoff-Witt basis), 
the Verma module Vo of the identity is found to have the 
(graded) character 



ch-V = - n - (1±5r " 3) 



(3)c 



(3.47) 



just as if all the modes L— n , G_„ + i/2 for 77 > (anti-) 
commuted (although they don't). Kac studied singular 
("null") vectors in such modules for the Virasoro and 
M = 1 superconformal algebras, and found in particular 
that the vacuum module, for which the conformal weight 
is 0, has a singular vector at level 1 in the Virasoro sec- 
tor, and level 1/2 in the superconformal current sector, 
for any value of the central charge c. The quotient by 
these singular vectors, which generate a submodule in the 
Verma module, leaves a module with exactly the charac- 
ter ch q V, and it is irreducible for generic values of c. 
That is, L—i and G-1/2 annihilate the vacuum in the 
quotient module. (For the Virasoro part, the fact that 
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the vacuum is annihilated by L_i means it is transla- 
tion invariant in the radial quantization point of view.) 
This conclusion is precisely what we expected according 
to Theorem 1. 

To avoid confusion, we should emphasize that the VOA 
generated by a(z) and a(z)*, which describes the chiral 
algebra of the "edge theory" of our states (or would, if 
the bulk were gapped) is not superconformal; it does not 
contain a superconformal current of weight 3/2, as one 
may see from the character above in the neutral sector 
(like all other half-odd-integer weights, it is removed from 
the character ch V by the projection to the even sector). 
Of course, the original superconformal current G still oc- 
curs, but only as part of a charged field a. For the basic 
case v — 1/3, this operator has weight 3 and is bosonic. 
The chiral algebra does contain a copy of the Virasoro 
algebra. 

We can now draw a conclusion about the space I of 
zero-energy states of the special Hamiltonian. We know 
that it must contain the space of amplitudes In of the 
VOA as we have discussed. But the count of states shows 
that in fact, for any non-zero value of c, these spaces co- 
incide at infinite N, in a suitable sense, that involves a 
shift in both the charge and degree that goes to infinity 
in the limit. This implies that the spaces In and In are 
the same in a certain limit of large N and large degree 
d. But we can now prove 

Theorem 2: In the example of this section, the spaces 
of polynomials In (the amplitudes) and In (the zero- 
energy states) are the same for any particle number TV. 
Sketch of proof: / is dual to A modulo the ideal gener- 
ated by the three-particle ideal I3 = I n ^3 , which tells 
us what combinations of £_ n 's to set to zero. In other 
words, it is defined as all those symmetric polynomials 
(in any number of variables) satisfying the vanishing con- 
ditions on three particles. Suppose the ideal T contained 
some further elements not in the ideal generated by I3, 
and that these first occur for N — N\ > 3. Equivalently, 
for N > Ni, In would require some further vanishing 
conditions not accounted for in In- For such N the di- 
mension dim//v,d would be smaller than dim/jv,d at some 
d. But for any N, the ideal In is translationally invari- 
ant, L—iXn Q Xv; its generators are local expressions 
in a and its derivatives, evaluated at z, and the corre- 
sponding relations hold for all z. In other words, the 
vanishing conditions involve only differences of Na of the 
coordinates, not the center of mass of these (the differ- 
ences are the relative or internal variables that appeared 
earlier). If for any N, dim In is smaller than dim/jy 
at some degree due to additional vanishing conditions, 
these affect only the internal variables, and so the effect 
propagates to all larger degrees because of the center of 
mass degree of freedom, which contributes (1 — q)^ 1 to 
the character. It also propagates to larger N, because in 
terms of I, any additional relations can be simply multi- 
plied by a(z). Consequently the dimensions would obey 
dim I N,d < dim70v,d (strict inequality) for all asymptot- 



ically large N and d. This contradicts the fact that the 
asymptotics of the two are known to be the same. Hence 
the ideal X is generated by I3, and the spaces of ampli- 
tudes In coincide with the spaces of zero-energy states 
In for all N. The characters of the components of the 
principal subspace are ch g Wn = ch g In, which we have 
calculated. (For Wn, we define the character unconven- 
tionally, as ch g Wn = tr WN q L °~ Nh .) 

We want to emphasize that the technique used to prove 
Theorem 2 is very general. It can be applied to other 
examples in which some translation-invariant relations 
involving some number of a's and their derivatives are 
known, or equivalently given a special Hamiltonian, if 
two conditions are met. These are that we can count the 
number of zero-energy states in each degree d and par- 
ticle number N, and that these numbers agree with the 
vacuum module V of the CFT (or VOA) asymptotically 
in the standard limit for dealing with the edge states. 
When these conditions are met, the equality then extends 
to finite N. We point out that it has been a perennial 
problem when using the MR construction to show that 
the zero-energy states of some special Hamiltonian are 
precisely the spaces of amplitudes, and no larger. For 
the case of the RR states, the proof was constructed ear- 
lier [40j , but to the present authors it appears technically 
difficult, and it may not be possible to extend it to other 
cases. By contrast, the argument given here deals with 
the MR and RR states using only the counts of states 
that have long been known. In FJM, the same result is 
reached in their examples, but the last part of the argu- 
ment is not explicitly given. Our approach also differs in 
that we have methods to find all the zero-energy states 
explicitly, and not only to give an upper bound. These 
methods were completed for the RR states by the authors 
of Ref. [H, but they quoted the old result of Ref. HJ to 
obtain the equality of the spaces of polynomials, rather 
than obtaining that as a byproduct of the construction. 

Next we turn to the inverse problem, finding special 
Hamiltonians such that the zero-energy subspaces are 
some desired spaces of functions, such as the amplitudes 
In or Jn of a given MR construction. This is essen- 
tially the problem of finding "vanishing conditions" that 
characterize these spaces. If the spaces of amplitudes 
In or Jn can be defined as symmetric polynomials that 
satisfy some "vanishing conditions" as ZiS come to the 
same point, we may construct a Hamiltonian that an- 
nihilates these functions and gives positive eigenvalues 
to functions not in the subspace, by forming projection 
operators that project onto functions (or states in the 
quantum mechanical Hilbert space, of lowest LL states 
of N particles) that are not in the subspace. Now the 
amplitudes form a subspace, and functions linearly in- 
dependent from them lie in the quotient space by that 
subspace. There is no natural subspace of representa- 
tives for this quotient space, in the absence of an inner 
product. In the quantum mechanical situation, we may 
use the L 2 inner product on lowest LL states of the plane, 
sphere, or torus, depending on our own interest. Hence 
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the form of the projection operator, onto some functions 
(in N' — 2, 3, ...variables) linearly independent from 
the "allowed" subspace, necessarily depends on the sys- 
tem geometry, and even on its size (Nj,). Hamiltonians 
that are a sum of such terms (necessarily invariant under 
permutations) with positive coefficients (the magnitude 
of the coefficients is not important for our discussion) 
have been in use in the QH effect for some time (see e.g. 
d, 03, [IH IHj])- Because the relations to be imposed on 
the functions are local, the projection operators in these 
special Hamiltonians are short range, even if they involve 
interactions among more than two particles. 

In the known examples, the interactions represented 
by the projection operators involve a bounded number 
of particles, and bounded degree for the relations (total 
number of derivatives appearing in the vanishing condi- 
tions) for each number of particles. If these conditions 
are not met, then the special Hamiltonian would not be 
local. (Higher degrees, which must necessarily occur if 
more particles are involved, imply a more extended in- 
teraction). Actually, this need not be strictly true. We 
only require the projection operators to have positive co- 
efficients, while the ground and quasihole states would 
be of zero energy. The higher degree, and so longer 
range, terms could be given numerically smaller coeffi- 
cients, and then the Hamiltonian can be relatively short 
range (or effectively involve power laws in distance) . But 
it seems preferable in aiming for physical examples that 
the boundedness in particle number and degree of the 
terms be respected as far as possible. 

Because of the duality discussed in Sec. MI B[ we can 
describe a set of spaces of polynomials (one space for each 
N) in terms of the algebra A and and a translationally- 
invariant ideal X, with the action of the positive cur- 
rents modes (and hence, symmetric polynomials) and the 
translation generator as before, even in the absence of a 
MR construction from a CFT. (It is less clear if we should 
also insist on an action of the positive half of the Virasoro 
algebra.) Then the space of polynomials that satisfy the 
vanishing conditions obeys the duality I = A/X, This is 
the point of view we will take for a while. 

In terms of the discussion of the spaces /, and the 
relations X in A, this means that the former are deter- 
mined by a set of translation-invariant vanishing condi- 
tions, that should preferably be bounded in both Na and 
d, and consequently finite. In terms of the relations in 
the ideal X, this property admits a clear mathematical 
definition: there is a finite set of elements i\ (indexed 
by I), each ii lying in a subspace Xd lt Nn sucn that X is 
generated from by taking linear combinations with 
coefficients in A and using the action of positive current 
modes, and the translation generator L-\. (We should 
point out that L_i acts on £_„ in A in the same way as 
on a_„, except for the shift in index by h, and so raises 
the total degree of excitation J2i n i-> as it is the transpose 
of the derivative operator in the dual space of symmetric 
polynomials.) It is far from clear if this finiteness prop- 
erty will generally be the case for the VOAs of interest 



here. This question now becomes an important issue to 
be studied in future work. 

To describe the structure of the special Hamiltonians 
in practice, some additional notation will be useful. To 
present the information about what terms appear in the 
special Hamiltonian, we can form a character for the part 
involving interactions among N a particles, in which the 
coefficients are the dimensions of the space onto which 
the projection operators project in each degree. For ex- 
ample, for the special Hamiltonian for the generic SCFT 
model, the series is 

1 + q 2 + q 3 + q 4 + q 5 + q 6 (3.48) 

for the three-body terms, and zero for Na ^ 3. This se- 
ries is the character of the space I3 , omitting the factor 
(1 — q)^ 1 which describes the uninteresting and ubiq- 
uitous center of mass degree of freedom that relates to 
translation invariance. Comparing with the series in Eq. 
(| 1 . 13[) for the internal part for the space of all symmetric 
polynomials, we see that the first term allowed is one of 
those in degree six, followed by all those of higher degree. 
As these series are indeed finite at least in the examples 
we will discuss, it is simpler to present this information 
than the character for the "allowed" functions in I. (Re- 
call that ch g / 3 = ch 9 A 3 — ch q X 3 .) Thus this type of 
series allows a succinct description of the terms in the 
Hamiltonian, though without giving the precise vector 
onto which projection is performed, which is important 
for the complete description when there is more than one, 
as occurs here for degree six for three particles. We seek 
to identify a set of generators (in the sense above) for X, 
and the use of a corresponding term in the Hamiltonian 
for each one should produce as the zero-energy subspace 
I N for all N. 

Because we consider finite sets of generators, it makes 
sense to add generators one by one, that is to add addi- 
tional translation-invariant terms to the special Hamilto- 
nian. In addition we must preserve the property that X 
is a module over the (dual action of the) symmetric poly- 
nomials or positive current modes, if we are to have any 
hope of finding polynomials that are zero-energy states 
for all N, or of agreeing with some MR construction from 
CFT. On X, the current modes j n act by lowering the ex- 
citation degree, as they are the transposes of the sums of 
powers s n that raise the degree of the symmetric polyno- 
mials. 

The special Hamiltonian for the SCFT model above, 
like other known simple examples, satisfies these proper- 
ties in a basically trivial way. Further cases will require 
less trivial checks that the ideal is a module over the pos- 
itive current modes. To find examples, we return to the 
observation that the space of amplitudes Jjv , constructed 
from inner products, which correspond to the entangle- 
ment subspaces of the (very large) ground state, may be 
smaller than the space of amplitudes ijv constructed from 
the dual space, when the latter contains singular vectors. 
The "missing" vectors in Jn give us an opportunity to 
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add elements to I, so that the missing vectors are those 
forbidden by the additional vanishing conditions. 
For the Af = 1 SCFT, it is known 50] for 



3 / 2(p-p'f 

c = ~ 1 ; 

2 V pp' 



(3.49) 



where p, p' are integers, p— p' is even, and p and {p—p') /2 
are coprime, that there are additional singular vectors in 
the modules, in particular in the vacuum module. This 
implies that they will occur at some corresponding de- 
gree in the spaces In also, at least for sufficiently large 
N. Setting all singular vectors to zero results in the su- 
perconformal minimal models, denoted SM(p,p'). One 
case is the tricritical Ising model example, which how- 
ever still seems to be difficult. Here we will pursue a 
simpler example, studied in FJM (though not from the 
SCFT point of view). Like the tricritical Ising model, 
the M(3,8) Virasoro minimal model (which has central 
charge c = —21/4) contains the J\f = 1 superconformal 
current at a corner in the Kac table, (2, 1) in the M(3, 8) 
case, and hence M(3,8) = SM(2,8). FJM showed that 
in this case, the ideal for three particles has the character 

(g)lch g X 3 = l + q 2 + q 3 +q 4 + q 5 +q 6 + q 8 , (3.50) 

and that I3 generates I. The corresponding vanishing 
conditions thus determine spaces In — In of polynomi- 
als for each N which are those of the VOA of the SCFT 
model, with the singular vectors set to zero. Compared 
with the special Hamiltonian of the generic superconfor- 
mal model, an additional translation-invariant term has 
been added in degree 8, for three particles, and nothing 
else. In the next section, we analyze this model, and 
some others, in detail. These models provide non-trivial 
examples of the structure we have just discussed. 



IV. ANALYSIS OF SPECIAL HAMILTONIANS 
FOR M(3,p) SERIES 

In FJM [23j], the example described at the end of the 
previous section was analyzed in detail. It is just one in 
a series in which the CFT used is the M(3,p) minimal 
model, p = 4, 5, 7, . . . (p is not divisible by 3), with the 
central charge of the M(3,p) minimal model 



c= 1 



2QP-3) 2 
P 



(4.1) 



The field ip is the (2, 1) field in the Kac table, with con- 
formal weight hjf, = (p — 2)/4. We concentrate on the 
basic cases of bosons, for which the filling factor is again 
v = l/(2/i^,). These examples of the MR construction 
have also been described as (k,r) — (2,p — 2) models; 
here as before k — 2 means they are paired states, while 
r = p— 2 is the degree of the lowest non- vanishing allowed 
behavior for k + 1 = 2 variables. 

The analysis of FJM begins with the ground state 
wavefunction (symmetric polynomial) they call ip^. They 



show that the space In is generated over the symmet- 
ric polynomials by this and £i(f3), and compute the 
character, as presented above for M(3,8). Taking the 
three-particle ideal I3 as determining vanishing condi- 
tions, they then analyze the zero-energy subspaces for all 
N. Because a bound on the characters of these agrees 
asymptotically with a so-called fermionic character for- 
mula derived from that for the M(3,p) minimal model, 
they conclude that I3 generates X, and In — In for all 
N. They also calculate the characters for functions on 
the sphere, corresponding to counting quasihole states. 

We will not repeat more details of FJM's analysis here. 
That is because the remaining methods do not obviously 
generalize. We do two things: we construct explicit sets 
of wavefunctions for all the zero-energy states, along sim- 
ilar lines as above for the generic superconformal model; 
and we describe the three-body interaction Hamiltonian 
for which these are the zero-energy states, for both the 
plane and the sphere. For completeness, we do the same 
for the other members of the M(3,p) sequence. 



A. M(3, 8) wavefunctions in the plane 



To find the wavefunctions, we make use of FJM's re- 
sults. The vanishing conditions on the TV- variable func- 
tions are incorporated by requiring that any function oc- 
curring in the particle-partition decomposition for any 
Na = 3 variables must lie in the span of <p$ and £1(^3) 
times symmetric polynomials, that is in I3 = I3. They 
use a filtration obtained as in Sec. [IT] by setting variables 
equal in pairs. A key point is that in the present case, 
because of the additional vanishing condition, the ker- 
nel of this map applied to three particles lies not only in 
DnIn R , it lies in its subspace DnI^j where I N { is the 
space of symmetric polynomials that vanish as degree 
three or higher when any three variables coincide (the 
ground state or lowest degree polynomial in this space 
has been called the GafFnian function (24[). These are 
related to the M(3, 5) minimal model, and are again ex- 
plicitly known [24J. 

Combining this observation with the construction of 
functions as before, we find a complete and linearly- 

— fSl 

independent set of polynomials in I N , essentially by re- 
placing elements of the ideal /™ 2 , that was used be- 
fore, with elements of iN-irlt ■ Explicitly, with notation 
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as before, we replace (4.9) with 



B. M(3, 8) wavefunctions on the sphere 



l(i,i)>(i,iO 



>< n n n^ws^s?)- 

a'=2,3 (l,i)>(o',i') J=l 

>< n n^^ ( j +1 ;*a } )- 

(2,t)>(a',i') j=l 

x n (4?-^) 2 - 

(3,t)>(3,t') 



a i=l 



(a) _ 



Here the (a, i)th cluster center of mass is Z 

— y^ r " , zj-"' and the only difference from the generic 
SCFT case is in the third line, where we define 

Xz\zi,z 2l z 3 ) = (z 1 - z 3 ) 2 (z 2 - z 3 ), (4.3) 
X Gi ( z i, z 2;z 3 ,z 4 ) = x° f (zi, z 2 ; z 3 )x 3 f (z 2 , Z\\ z 4 ). (4.4) 

Note that this only affects couplings between pairs of 
particles of type a — 2, and between a pair of type two 
and a particle of type three. 

The degree of each function, omitting the elementary 
symmetric polynomial factors, is now 

3 1 

-N(N - 2) + m\ + m 2 m 3 + -mj + m 2 + m 3 . (4.5) 
Hence the character is identical to that given in FJM: 



E 



m 2 ,m3>0:2m2+m3< N. 

(-1)'"3 =( _1)« 



(q) »-2m r m 8 (q) m , (q) 

7 



(4.6) 



For N — > oo we can add the odd and even sectors to 
recover the SM (2, 8) vacuum character in fermionic form 

da 

MoocMtf? , (2 < 8) = (9)., U\V^ + q^ q V ( f^ 



E 



m 2 i m 3^0 



(q)m 2 (q)r 



the first few terms of which are 
7 (i,i) 



(q)ocCh q V ( f[*' 8) =l + q 2 +q 3 + 2q 4 + 2q 5 + Aq 6 



(4.7) 



(4.8) 



(<z)ooCh 9 V/ ( f^ 8) = q 3 / 2 (l + q + q 2 + 2q 3 + 3<? 4 + 4q 5 + . . . 

(4.9) 

These should be compared with those for the vacuum 
Kac module of the generic case. 



As for the generic superconformal model, the wave- 
functions on the sphere can be found with only slightly 
more effort. The degree in each coordinate variable to be 
N<f, — (6N + n)/2 — 6, as this is determined only by the 
"leading order" / "ground state" clustering behavior pa- 
rameters (k — 2,r = 6). Repeating the power counting 
in the residues of the wavefunctions (|4.2p . we find that 
the required number of quasiholes of each type are 



n\ = n, 

n 2 = n — 2m 2 — m 3 — 2, 
n 

n 3 = - - m 2 - m 3 - 2. 



(4.10) 
(4.11) 

(4.12) 



The factors involving particles of type one are unchanged; 
we only need to homogenize the coupling (|4. 3[) between 
a half-broken pair of type a = 2 and a particle of type 
three: 

xf f (zi,z 2 ;z 3 ;wi) = (zi~z 3 ) 2 (z 2 ~-z 3 )(z 1 ~w 1 )(z 2 ~w 1 ) 2 . 

(4.13) 

Again, because n 2 = 2n 3 + m 3 + 2, there are always 
enough quasiholes of type two to do this. The basis func- 
tions on the sphere are therefore 



4 n x(^S;4;^K) 

l(l,»)>(o,i') 



X3\\l , z i,2, Z i',V W 2i'-l, W 2i' ) 

(l,t)>(3,t') 

x n *S) 

(2,i)>(2,i') 

x n 

(2,»)>(3,i') 
(3,i)>(3,i') 

xnnn^^f-^)}- ^ 

a i,j I ) 

The ranges of the quasihole factors in the last line are 



= 2m 3 + 1, 

= 2m 3 + n 3 + 3, 

= m 3 + 1, 

= (n 2 + m 3 + 2)/2, 

= 1. 



,2m 3 + n 3 + 2, 
,111, 

,(n 2 + m 3 )/2, 

,"-2, 
,"-3, 



(a = l,j = 1), 
(a = l,j = 2), 
(a = 2,j = l), 
(a = 2,j = 2), 
(a = 3, j = 1). 
(4.15) 



It may be verified that the character on the sphere is 
still of the form 



tr q 



n 



m a + n a 



(4.16) 
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with the n a now given by (|4.10p - (|4.12[) ; this is equivalent 
to theorem 5.14 in FJM (where the variables + 1, N 
are denoted N, n). Taking the planar limit n — > oo with 
N fixed, we recover (|4.6p . Furthermore, the quasihole 
multiplicity for fixed positions is now 



n — mi — m 3 
m 2 



n/2 — vri2 
m 3 



E 

m2,m3>0: 
N-2m 2 -m 3 >0, 

(-l)-3=(-l)«. 

(4-17) 

the domain of which is clearly bounded and independent 
of N for N > n (hence the sum itself also has these prop- 
erties). Even though this multiplicity is now bounded, 
unlike in the generic superconformal case, we still expect 
this model to be gapless in the thermodynamic limit be- 
cause it is a non-unitary CFT. 



C. M(3,p) series: general p 

As found by FJM, the structure discussed above for 
the M(3, 8) case holds for all M(3,p) models with p > 5. 
That is, the kernel of each map C\ for 1$ is D^I^~ 3 ^ 

(and note that f§ = lj? R ). Then the spaces of zero- 
energy states for the M(3,p) models can be found recur- 
sively for all p not divisible by 3. Defining 



(4.18) 



there will be a = 1, . . . , s types of particles with r a = 2, 
representing pairs in various stages of "brokenness," and 
one type of unpaired particles with a — s + 1, r a = I. 
Schematically, the plane wavefunctions take the form 



Mil 



A \ z iS ) z i,2 i *z',H z i' .2) 

(a,i)>(a,i') 



(a)f (a) (a) _ (ct')\ 

*->A A 1 A I 1 1 &AI 



(a,i)>(a',i') i=l 



X D 



_ /f (/^=a+l,..., S +l)l\ 

«nn^({^})" 



(4.19) 



a ;=i 



All pairs of rows of type 1 < a < s are coupled by a 
function 

x {a) ( z i, z 2;z 3 ,z 4 ) = 

E 4" ) (( 2:i_2;3 )( Z2 ~ 2;4 )) P 3 " +1 J (( Z 1 _Z 4)(^2-^3)) J , 



where the Cj are constants which are partially fixed by 
the requirement that the symmetrization of this function 
reproduce the four-point conformal block: 

S*X (a) (zi,z 2 ;z 3 ,Z4) 

= (a( Zl ) ■ ■ ■ a(z4)) M(3p _ 3(a _ 1))m{1) . (4.21) 

Because these are Virasoro minimal models, the four- 
point block may be evaluated in terms of the ordinary 
hypergeometric function as [23] 



S z X^ a \zi,Z2;Z3,Zi) = ((zi - z 3 )(z 2 - Z4)) 1 



-3q+1 



x . 2ri {a-^,3a-p-l;2a-^-,x) . (4,22) 



where the cross-ratio x = (zi — z±)(z2 — z 3 )/(zi — z 3 )(z2 — 
Z4). The leading-order behavior in z\,Z2,z% as Z4 — > 00 

then defines X3 ( z i, z 2', 23). 

We use the term "schematic" in referring to the wave- 
functions (|4.19l) because x 1S n °t uniquely determined 
by Eq. (14. 21)) : there exist non-symmetric functions of 
the form (|4.20l) which vanish under symmetrization but 
whose contributions do not vanish when used to con- 
struct a wavefunction of the form (|4.19[) with N > 4 
particles. The ambiguity can of course be resolved by an 
examination of the singularity structure of the four-point 
block (■ij'(zi) ■ ■ ■ 0(z4))jw(3,p-3(a-i))j as was done for G in 
Ref. [22( , but this must be done on a case- by-case basis as 
the conformal weight of ip (and hence the number of sin- 
gular terms in the t/j ■ t/j OPE) increases with p; therefore 
this is the most we can say about the form of the wave- 
functions without restricting ourselves to specific cases. 

For general p, characters on the plane and sphere may 
still be computed without explicit basis functions; this 
was done in section 5 of FJM. We restate their results in 
a form more conventional in the FQHE literature [iil. I52I] : 
let 



(4.23) 



denote the quasihole degeneracy at fixed positions in the 
M(3, p)®J7(l) theory on the sphere with n quasiholes and 
F not-fully-paired particles (in the notation of (|4.19l) . 



F = J2 a= 2 r a m a)- The relationship between N,n and 
Nj, is 



N 



[p] _P 



-N - 



n 
2 



(P-2). 



(4.24) 



FJM's equations (5.5), (5.6) then give a recursive expres- 
sion for the quasihole degeneracy, analogous to the result 
obtained in Ref. [44[ for parafermions: 



(4.20) 



E 



F'=F 



3d 2 



F-F' 



(4.25) 



p-3 
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with n' = n — F — 2 and initial condition 
_ff + (p-4)(f-l) 



F 



for p = 4, 5. (4.26) 



FJM solve this recursion in their theorem 5.14; in partic- 
ular the numbers of quasiholes at each stage are 



n\P ] = n 



2(a- 1) - ^(min(a,/3) - l)r p mp (4.27) 
0=1 



for 1 < a < s and 



[ P ] n 



^ s+l 



0=1 

+ i((pmod3)-2)(m s+ i-2). (4.28) 



D. Explicit special Hamiltonians on the sphere and 
the plane 

So far we specified the allowed behavior for any three 
particles as that belonging to the spaces Jg for three par- 
ticles. As these functions are explicitly known, we have 
only to construct translation-invariant projection opera- 
tors onto the orthogonal complement of these spaces, and 
the sum of these any positive coefficients gives a special 
Hamiltonian that produces the allowed functions as zero- 
energy states. Unlike the allowed functions, the forbid- 
den behavior, which must be orthogonal to the allowed, 
depends on the geometry of the system through the inner 
product in use. We give results for the sphere, for which 
the measure in stereographic coordinates (by mapping to 
the plane) is 



n(z) d 2 z 



dz d~z 



(1 + |^| 2 /4) 2 ( 1+JV */ 2 )' 



The results for the plane may be obtained by keeping 
only the leading order terms as N$ — > oo. 



FJM's results determine the character of the ideal X, 



[p] 



for the M(3,p) case as 



(q)ich q X 3 



[p] _ C 1 



(l-q 2 )(l-g3) 



(4.30) 



The series on the right-hand side terminate with a finite 
number of terms. The coefficients of these terms tell us 
how many orthogonal states must be forbidden in the 
polynomial behavior for each p. 

As a basis, we use the translationally-invariant sym- 
metric polynomials {e„}, 2 > n > N, defined in Ref. [3(| 
in terms of the elementary symmetric polynomials e„ as 



e« ({z»}) = e„ < z 



1 



N 

nY, z o 



(4.31) 



We specify the orthogonal allowed and the forbidden be- 
havior for each degree d in which the projection operator 
has non-zero image and kernel among the functions of 
that degree, since otherwise an arbitrary basis for that 
degree can be used (or omitted entirely). Recall that the 
lowest allowed degree is r = p — 2. The functions have 
not been normalized. The first examples are as follows: 
p = 7 (r = 5): At d = 6, the allowed function is 
9e§ — 2e§; the forbidden (to be projected onto in the 
special Hamiltonian) is 

27 (9iV| - A8N 2 + 73A^ - 40) e§ 

+ (63iV| - 471iV$ + 1156iV - 880) e|. (4.32) 

p = 8 (r = 6): At d = 6, allow Lpf oc 9ef, - e\\ forbid 

27 (9iV| - 607V| + 127iV - 100) e§ 

+ 2 (45iV^ - 4087V 2 , + 11997V - 1100) e\. (4.33) 

The behavior at d — 7 is allowed, and at d = 8 allow 

e^cpjf' oc 9e§e 2 — ef ; forbid 

54 (9N% - 69N 2 + 160A^ - 140) ^e 2 

+ (99iV^ - 9487V 2 , + 29997V - 3010) e\. (4.34) 

p = 10 (r = 8): At d = 8, allow tp [ ™ ] oc -18e|e 2 + ef ; 
forbid 



27 (9iV| - 93A^| + 316JV - 392) ^e 2 

+ (63iV| - 7867V, 2 . + 3217/V - 4214) e 2 . (4.35) 



At d= 9, allow (li - 16/3si)4 101 oc 3e| - 2e 3 e|; forbid 

27 (9JV$ - 84AT 2 , + 235iV - 216) 

+ (63 iV^ - 8047V 2 , + 33977V - 4392) e z e\ . (4.36) 

[10] 



(4.29) At d = 10, allow e 2 co> 1 oc -18epl + e|; forbid 



135 (9iV| - 102iV| + 3737V - 504) e%e% 
+ 2 (ll7iVj - 15427V 2 , + 67457V - 9576) e|. (4.37) 

Both behaviors at d = 11 are allowed, and at d = 12 
we have a three-dimensional space; we must allow the 
two-dimensional space spanned by the (nonorthogonal) 

behaviors e3(£i — 16/3si)(/53 10 ' oc 3ef — 2e|e| andef,^ 10 ' oc 
— 18e|e| + e|, and forbid the orthogonal complement 

3645 (81JV| - 17557V| + 14607iV| 

-601177V| + 1342847V 2 , - 169756iV + 110880) e\ 

+ 270 (1053/V^ 5 - 30834iV| + 369873iV^ - 2330376A^| 

+8190988iV| - 15458176/V + 12418560) e 2 ,^ 

+ 2 (20007/V^ 5 - 641979/V^ + 84799177^ - 58925481/V| 

+2268509567V, 2 . - 458O5366OiV + 378655200) e|. 

(4.38) 
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None of the coefficients above factors as a polynomial 
in N^, and using different bases for the polynomials [such 
as Jack polynomials at the appropriate parameter value 
a = — (k + l)/(r — 1)] do not yield noticeable patterns in 
the expressions involved. 



V. CONCLUSION 

To conclude, we have explored the relations among sev- 
eral spaces of wavefunctions that can be constructed for a 
system of trial wavefunctions, especially when those can 
be constructed from a CFT. These were: 1) the entan- 
glement subspaces of the ground state, for Nj± particles 
out of N; 2) spaces of amplitudes (conformal blocks) con- 
structed from a CFT; and 3) spaces of zero-energy states 
of a special or projection operator Hamiltonian. In fa- 
vorable cases, we can show that all of these consist of 
exactly the same wavefunctions. Even when they do not, 
it may be possible to add terms to the Hamiltonian and 
restore agreement, in terms of a "quotient" CFT. This 
is a finite-size version of a bulk-edge correspondence, in 
current parlance. 

We have not addressed the issue of the energy gap 
above the zero-energy states for the special Hamiltonian. 
It is expected that when the CFT is non-unitary or non- 
rational, the model will be gapless in the thermodynamic 
limit 14]. We want to remark here that our analy- 
sis makes it clear that there are many (many-particle- 
interaction) terms that can be added to such a Hamil- 
tonian, such that the zero-energy states are unchanged. 
The non-zero energy states would be affected, however. 
But it should be the case that no such terms can produce 
a gap in the bulk spectrum. This and the nature of the 
gapless systems remains a topic for future investigation. 
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Appendix A: Proof of Theorem 1 

First, we prove that a(z) and a(z)* generate the VOA. 
We take as given these operators at any locations z, and 
use translation invariance of ope's; we must show that 
the current and Virasoro algebras are generated by these 
operators. To begin, we take the ope a{z')a{z) as z' — > z 



[z fixed), and extract the leading term, which contains 
the operator e 2lv ^ z ^ ^ times the leading operator in the 
rptjj ope; the whole operator is located at any z. Then 
we fuse this with another a, and repeat the process. By 
definition of the properties of after taking the prod- 
uct of k a's, we obtain the operator e %k vi z )l V 77 (times the 
identity in the ?p theory). Doing the same with a*'s pro- 
duces e~ %k ^ z - . The ope of these two operators in the 
charge sector (the ip CFT) contains the current operator 
j(z), and the ope of two js produces the stress tensor 
T v (z). It follows that all descendants generated by cur- 
rents ("current excitations") of the leading operator in 
the ope a(z')a(z) as z' — > z are also in the algebra, and 
similarly for the leading operator when another a(z") ap- 
proaches this, and so on. 

Next, we consider again the ope of a with itself. In 
addition to the leading term referred to above, it also 
contains subleading terms, each of which is a sum of the 
leading and sub-leading terms in the iptjj ope, possibly 
with current excitations, with weights that sum to the 
correct value corresponding to the overall power of z. 
The subleading terms consisting of current excitations 
only can be subtracted off, as they are already in the al- 
gebra. At each level (extra factor of 2), that leaves an 
operator containing descendants in the algebra generated 
by tj), times current excitations. At the first level, there is 
just one term, the first subleading part of the iptp ope (if 
it is non-zero), times e 2lip ( '/ v", which is therefore also in 
the algebra. At the next level, the current descendants 
of this can be subtracted off as well, possibly leaving an- 
other operator with only descendants in the ip algebra. 
Proceeding in this way, we find that all descendants in 
the ^{^^{z) ope, times e 2l¥>( - ^ v/ ^ 7 , times current exci- 
tations, are in the algebra. This can be repeated also in 
the ope of these operators with a further a(z'), and so 
on. Then we obtain recursively all of the VOA V, de- 
fined in Eq. (|3.18j) . In particular, the stress tensor T^(z) 
is eventually generated. 

We can now prove Theorem 1. We observe that the 
iterated ope of k — 1 a's, taking the leading term in 
each ope, produces e 1 ^' 1 ^^^-^* (z), by definition be- 
cause a further ope with a gives e lkiplyZ ^ ^ . The ope 
with e- lk v\ z )l-J» gives a* as the leading term. Hence the 
VOA V is generated by a(z) and e^M 2 )/^ which is 
the second statement in Theorem 1. 

To obtain the first statement, we consider the opera- 
tors e l7V ^ 2 )/^, for N divisible by k, which (for z = 0) 
correspond to the states |0jv). Acting on one of these 
with a's gives the space V N , by definition. An ope with 
e -ifcip(z)/ % /l7 j owers jy k. Hence the union of the spaces 
of states (or corresponding operators) V~ N over all non- 
negative multiples of k, or the limit limj\r_ i . 00 V~ N using 
the embeddings, gives all of the VOA V. That is the first 
statement in Theorem 1. 
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